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Abstract 

We present a necessary and sufficient condition for three qutrit density matrices to be the one- 
particle reduced density matrices of a pure three-qutrit quantum state. The condition consists 
of seven classes of inequalities satisfied by the eigenvalues of these matrices. One of them is a 
generalization of a known inequality for the qubit case. Some of these inequalities are proved 
algebraically whereas the proof of the others uses the fact that a continuous function of the state 
must have a minimum. Construction of states satisfying these inequalities relies on a representation 
of the convex set of the allowed set of eigenvalues in terms of corner points. We also present a 
result for a more general quantum system concerning the nature of the boundary surface of the set 
of the allowed set of eigenvalues. 

PACS numbers: 03.67.Mn, 03.65.Ud, 02.40.Ft 
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I. INTRODUCTION 



Recently there has been much activity in clarifying the nature of entanglement in quantum 
systems of finite degrees of freedom in connection with quantum computing, and there have 
been some results about properties of reduced density matrices of a pure state. In particular, 
a necessary and sufficient condition has been found for n one-qubit density matrices to be 
the reduced density matrices of a pure n-qubit quantum state 0] for all n. [The phrase 
"reduced density matrix" is used to mean "one-particle reduced density matrix" . It will be 
abbreviated as RDM below.] In this paper, as a continuation of this line of investigation 
we find a necessary and sufficient condition for three one-qutrit density matrices to be the 
RDMs of a pure three-qutrit quantum state. (Throughout this paper "a state" means "a 
pure state" unless otherwise stated.) 

We consider the space of states spanned by the basis vectors \iik) = \i) ® \j) ® \k), 
where i, j and k are 1, 2 or 3. The one-particle states |2) and |3) are assumed to be 
orthonormal. A general state can be given by 

3 3 3 

i=i j=i k=i 

We require that (\E' | vp) = J2ijk kuAiP = 1- Thus, the state |\E') is a normalized three-qutrit 
quantum state. The three RDMs are given by p''^) = '}2j,k CijxCi'jK, pf) = 12i,k cijKCifW 
and p^^i!, = J cjjkCijk'- We denote the eigenvalues of the ath RDM, p[f , by A^"'', A2"^ and 
A^"^ with the ordering aS''^ < A^"^ < A^''^. We now state the main theorem of this paper. 

Theorem 1. A necessary and sufficient condition for nine nonnegative numbers X'f\ i = 



1, 2, 3, a = 1, 2, 3 satisfying A^ < A2 < A3 and A^ + Ag + A3 = 1 to be the eigenvalues 
of the three RDMs of a three-qutrit quantum state is given by the following inequalities: 

+ aS'^) < Af + Af) + A^^) + X['\ (1) 

A^'^) + aS"^ < A? + + + \['\ (2) 

A^'^) + A^") < Af + Af) + A^^) + X'i\ (3) 

2A^'^) + aS'^^ < 2A?) + + 2\f + \^\ (4) 

2Ai'^) + A^") < 2A^'^ + \f + 2Af^ + \^\ (5) 

2A^'^) + A^'^^ < 2Xf + aS') + 2A^^^ + A^^^ , (6) 

2A^'^) + A^'^^ < 2Af ^ + A^') + 2A^^^ + A^ ^ , (7) 



where {aba) are any permutations o/(123). 

If A^"'' = for all a, then by writing A3"'' = 1 — Ag'^'*, these inequalities can be expressed in 
terms of We find that the only nontrivial inequalities resulting from Theorem are 
X^^ < X^ + X'^'^ for all permutations (aba) of (123). This conclusion agrees with Refs. [H, 0|- 
In Ref. j2[ Bravyi solved the corresponding problem for the system ® ® C^. His result 
can be used to find a necessary and sufficient condition for three density matrices to be 
RDMs of the quantum system (g> C^. This condition can be compared with that 
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obtained from Theorem^ by letting = = 0. These conditions can be shown to be 
identical. 

The rest of the paper is organized as follows. In Sec. Owe represent the set of the allowed 
set of eigenvalues as a convex set determined uniquely by the corner points. In Sec. IIIII we 
prove a generalization of inequality (Q). Inequality itself follows from this as a corollary. 
In SecHVlwe provide a proof of inequalities (j21) and 0. Our proof of inequalities (H))-® is 
algebraic. The proof of the other inequalities uses the fact that a continuous function of the 
(ordered) set of coefficients (cijk) with the constraint ^jj^ iQjfcp = 1 must have a minimum 
because the domain can be identified with the 53-dimensional sphere of radius one, which 
is compact. We explain our method in Sec.|V|and use it to prove inequalities (jlj) and (0) in 
Sec. IVII In Sec. I VI II we prove inequalities (jH)) and ((Tj). [Our proof of inequality (jH)) is rather 
lengthy.] We then construct states for each set of eigenvalues satisfying (H}-© in Sees. IVIIII 
and IIXI The corner-point representation given in Sec. |n] is useful in this construction. We 
conclude this paper in Sec. |X|by showing that the set of the allowed set of eigenvalues for a 
more general quantum system is likely to be bounded by hyperplanes. 

II. THE CORNER-POINT REPRESENTATION OF THE SET OF THE AL- 
LOWED SET OF EIGENVALUES 

Let us write a set of possible eigenvalues as 

[^1 \ ^3 ' ^1 ^2 ^3 > ^1 ^2 ^3 J • W 

Since the combinations 001 and 0|| and ||| will appear often, we abbreviate them as 
O = 001, A = O^l and i? = A set of eigenvalues can be regarded as a point in a 

six-dimensional space with coordinates A^"'' and a = 1,2,3. For this reason we call a 
set given by (jS} an eigenvalue point, or an E-point in short. In this representation the set of 
I5-points defined in Theorem ^ is a convex set bounded by hyperplanes. Therefore this set, 
which will be denoted by S, can also be specified by giving all the corner points, i.e. the 
points on the boundary of the set where there is no tangent line to the boundary. 

If we impose only the conditions < A^"^ < A^'*^ < A^''^ and A^''^ + A^"^ + A^''^ = 1, then 
the set of all allowed E-points are bounded by the hyperplanes A^"'' = 0, A^*^^ = Ag"'' and 
A2"'' = Ag"'' = 1 — a["'' — X'^\ a = 1,2, 3. The set of E-points satisfying these conditions is of 
the form L x L x L, where L is given by the triangle shown in Fig. 1. 
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The corner points of the set L x L x L are 



[0,0,0], [0,0, A], [0,0, B], [0,A,A], [0,A,B], [0,B,B] 
[A, A, A], [A,A,B],[A,B,B], [B,B,B] 



(9) 



and the E-points obtained from these by qutrit permutations. (For example, the E-points 
[A, B, A] and [B, A, A] are obtained from [A, A, B] by qutrit permutations.) 

If Ag"^ = 1 and A2"^ = A^"'' =0 for some a, the state reduces to a two-qutrit system. It is 
well known that the two RDMs must have the same eigenvalues for any two-particle state. 
It can readily be verified that this fact is compatible with Theorem ^ In fact inequalities 
and © imply this fact and the rest are satisfied if Ag"'' = 1, Xf^ = X\'^\ i = 1, 2, 3 for some 
(abc). This implies that the E-points [O, O, A], [O, O, B] and [O, A, B] are not in the set S. 
All the other E-points in the list Q satisfy the inequalities of this theorem. Therefore they 
will remain corner points of S. 

Let us give the corner points of 5*. This will be useful when we construct quantum states 
satisfying the inequalities of Theorem ^ in Sees. IVIIII and IIXI 

Theorem 2. The convex set S of the E-points characterized by Theorem^has the following 
corner points and those obtained by qutrit permutations from them: 



Since the E-points in this list that are corner points oi LxLxL, i.e. the first seven E-points, 
remain corner points of S, we only need to show that the corner points of S that are not 
corner points of L x L x L are given by the last four E-points in (|Tn|l and those obtained by 
qutrit permutations from them. Since each of these new corner points must satisfy one of 
equalities (P)-®, i.e. equalities obtained by replacing inequality signs by equal signs, the 
following proposition, which lists all the corner points on the hyperplane given by each of 
equalities (H))-©, will imply Theorem El 

Proposition 1. The intersection of the set S and the hyperplane bounding each of the 
inequalities listed in Theorem^ is a 5-simplex. These 5-simplices are characterized by the 
corner points given as follows [for {abc) = (123)y.' 



m ■ [0,0,0], [B,0,B], [B,B,0], [A,0,A], [A,A,0], [5, Qif, Off] ; 

: [0,0,0], [B,0,B], [A,0,A], [0,A,A], [111,011 A]; 

m ■■ [0,0,0] [B,0,B], [A,0,A], [0i|, 011,5], [A,0\l\\l], [AMi,^]; 

: [0,0,0], [B,0,B], [B,B,0], [A,0,A], [A,A,0], 

© : [0,0,0], [B,0,B], [B,B,0], [A,0,A], [5, Qif, Oif] , [B,m,A]; 

m : [0,0,0], [B,0,B], [A,0,A], [A,0\l\l^], [A,|||,|lf], 

Q : [A,A,B], [B,0,B], [A,0,A], [oi|, 0||, S] , [AOHao]' [^111^]. 



[0,0,0], [0,A,A], [0,B,B], [A, A, A], [A,A,B], [A,B,B], [B,B,B] 

L '^33' 33J' L '^44'4 4 2J' L ' '663J' L '6 6 3'66 3J' 



(10) 
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Proof. Each set of E-points in this hst clearly forms a 5-simplex since there are six E-points 
each. It is straightforward to check that each simplex is nondegenerate, i.e. five-dimensional. 
To verify that these simplices form the boundaries, i.e. to show that there are no other E- 
points contained in the set 5" on each hyperplane, it is sufficient to show that the six sides of 
each 5-simplex are on boundary hyperplanes of S. This will be true if any five E-points chosen 
from the set of corner points of each 5-simplex are on a single boundary hyperplane of S. 
This fact can readily be verified. We list the hyperplanes bounding each 5-simplex together 
with the El-point to be removed. For example, for equality ((H), "[0,0,0] : X^2^ = Ag^^" 
below means that the E-points [B,0,B], [B,B,0], [A,0,A], [A,A,0] and [E,0||,0||] 

are on the hyperplane Ag^"* = X^^K 

Equality (0), (abc) = (123) 

[0,0,0] : A?) = Ai'); [B,0,B] : Af ^ = 0; 

[B,B,0] : aS^^ = 0; [A,0,A] : ©, (a6c) = (132); 

[A,A,0] : ©, (afec) = (123); [5,0i|,0||] : (gD, (a6c) = (123); 



Equality ©, (abc) = (123) 

[0,0,0] : Af =Af; [B,0,B] : A? ^ = 0; 

[A,0,A] : Af) = A^'^ [0,A,A] : ®, (a6c) = (321); 

[1||,|||,A] : Af) = 0; [iiiOif,^] : ®, (a6c) = (321); 

Equality ©, (abc) = (123) 

[0,0,0] : ^,{abc) = (123); [B,0,B] : aS'^ = 0; 
[A,0,A] : Af) = Af; [oi|,0||,5] : ®, (ate) = (123); 
[AOH,iii] ■■ ©,(afec) = (321); : Af ^ = 0; 

Equality (gj), (abc) = (123) 

[0,0,0] : A^') = A^'); [B,0,B] : @, (ate) = (231); 
[B,B,0] : (0),(ate) = (321); [A,0,A] : Af^ = Af ; 
[A,AO] : Af)=Af; : ©, (ate) = (123); 

Equality ®, (ate) = (123) 

[0,0,0] : A(^)=A«; [B,0,B] : Af^ = 0; 

[5,5,0] : ®,(ate) = (321); [A, O, A] : Aj'^ = A^'^ 

[i?,Oi|,0||] : X? = Xf; [B,l'^lA] : ©, (ate) = (123); 



Equality ©, (ate) = (123) 

[0,0,0] : A(^) = A«; [5,0,5] : aI^^ = 0; 

[AO, A] : Af) = Ar; [AOiiiii] : = 

[A, III III] : (ED, (ate) = (123); [A, Uf, 5] : @, (ate) = (321); 
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Equality ©, (abc) = (123) 

[A,A,B] : ^,{abc) = (123); [B,0,B] : A^^^ =0; 
[A,0,A] : Af) = Ar; [Oi|,0||,i?] : A?) = A«; 

[A OH, Hi] : Ar = Af; [A Mi 5] ^Af)=0- 

□ 

III. PROOF OF A GENERALIZATION OF INEQUALITY (HD 

Now we start proving the inequalities in Theorem ^ thus establishing that they are 
necessary conditions for the A,-"^ to be the eigenvalues of the RDMs of a three-qutrit quantum 
state. Inequality is a corollary of a more general result stated in Ref. 0. We present a 
proof of the latter. 

Theorem 3. Consider an n-qudit state given by 

l<ii,...,i„<d 

where the set , \2) , . . . , \d)} is an orthonormal basis for a d-dimensional vector space. Let 
, i = 1,2, . . . ,d, be the eigenvalues of the ath RDM with X^""^ < A^"'' if i < j . Then 

n-l d-l 

^^^(.)>^(n) + ^(n)^...^^W . (11) 

a=l 1=1 



Proof. The basis states can be chosen in such a way that the RDMs are all diagonal 
We work in this basis. Define 

n-l d-l 



a=l 1=1 

Note that 



\^ii---ia-iiia+i---in I • 

Let us write 

An) 



i«2---*n-ii ~ y ^hi 

If A}"'' > for all i, we have 



«l«2---«n-l 



^ ] '2ili2---in-li '^ni2-- in-l«' if ^ 7^ ^ . (l^) 



«l«2---«n-l 
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The latter follows from the fact that the nth RDM is diagonal. If A^""* = for some j, then 
the coefficients ajjj2...j^_ij are undefined. However, we can choose them so that Eqs. ()12p 
and p3p are satisfied for all i and i' even in this case. For a 7^ n we have 



d 

~ / ■^in / , \'^^\■■■^a-\^^aJr^-■■■1■n\ 

2n = l ,*a-l,«a + l,--- ,«n-l 



Then we find 



n-l d-\ 

1=1 i=l jl,i2---jn-l,«n 

where g{i\^ 12, . . . ,in) > 1 unless ii = i2 = ■ ■ ■ = in-i = d, and where g{d, d, . . . ,d,in) = 0. 
Thus 

d 



1 |2 






\^ili2---in\ 1 


— a 


dd---din 1 



= 1 - E as:v<.....<..„p, (14) 

in = l 

where we have used (fT^ . Eqs. (fT^ and (fT3|) allow us to consider aj^jj-'-in-ii components 
of orthonormal vectors labelled by i in c?""^ dimensions. This implies, in particular, that 
[cidd-di) add---d2, • • • , ddd - dd) IS part of a c?"~^ X c?""^ unitary matrix. Hence we have 

d 

J]|a,,...,.„P<l. (15) 

in = l 

The quantity X]l=i takes the maximum value under the constraint (fT^ if 

\o,dd---dd\ = 1 and add---di 
Then from ()14|) we find 



0'dd---dd\ = 1 and add---din = for i„ 7^ since A^"^ is the largest among A,^"^ by definition. 



G>l-Ai") = AS") + A^") + ---At\ 
as required. □ 

Inequality follows from this theorem by letting d = 3 and n = 3. 



IV. PROOF OF INEQUALITIES (HD AND © 

We work in the basis for which the RDMs are all diagonal unless otherwise stated. We 
let (abc) = (231) in inequality Q and (abc) = (132) in inequality © to find 

Af + Af < Af ) + Af + A« + A^^) , 
A« + A« < Af ) + Af + Af + Af ) . 
These are equivalent to the following proposition. 
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Proposition 2. A^^ - A^^ < A^^^ + A^^^ and Xf^ - A^^ < A^^^ + Xf\ 
Proof. First we note that 

Xf) _ Af ) < |ci23p + |C133P + |C223P + |c233p + |c323p + ^3331' • 

Define 

X = |ci23p + |Ci33p + |C223P + |c233p + |c323p + |c333p • 

Then it is sufficient to show that 

X < A« + A(^) , (16) 
X < X? + Af . (17) 



To prove (fTH|) we define bijk by Cijk = y A^^"*^^ bijk if A^^^'' 7^ for all i. Then we have 

3 3 

j=l k=l 

If xf'^ = for some i, then are arbitrary, but we can still choose them to satisfy these 
equations. Hence the following matrix is part of a 9 x 9 unitary matrix: 

— ( ^133 ^233 ^333 
\ &123 ^223 ^323 



Thus, we have the following inequalities: 

|&133r + |&233r+ |&333r < 1, (18) 

|&123r + |&223|'+|&323|' < 1 (19) 

and 

\h^?.? + \hi2^?<l foralH. (20) 

Note that 

X = X^\\hr,,\^ + |6i23n + A^'^(|6233r + \h22,?) + A^'^(|&333P + |&323n • (21) 

From (HH) and ^ we find 

I&I33I' + I&I23I' < 2 - (I6233I' + I&223I' + I&333I' + |&323n • 

By using this in ()2ip we obtain 

X < 2X? + (A^') - Al'))(|6233r + \b223\') + (A^^^ - A«)(|6333r + \h23\') 

< X^'> + A« , (22) 

where we have used ()20|) with i = 2 and 3. Thus we have established inequality ()16|). 

To establish ()17|) we first make a unitary transformation on the first and second rows of 
Ml and obtain 

^ I abi33 + Pbi23 a&233 + /3&223 a&333 + Ph23 
^ -/56l33 + a&123 -/?&233 + "^223 -/5&333 + a&323 
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where + = 1. We choose a and /3 so that —/5 6333 + 0:6323 = 0. If 6333 = C333/ y Ag^'* 7^ 
or 6323 = C323 / \j y^l^ 7^ 0, then a and /? can be chosen as follows: 

a = ' , (23) 

V IC333P + IC323P 

P = ^ . (24) 

V IC333P + |c323r 

If C333 = C323 = 0, then {a, (3) can be chosen, for example, to be (1,0). Then we have from 
the second row of Mi 

I - Pbi33 + a6i23p + I - /36233 + a&223p < 1 • 

By multiplying this by X^2^ and remembering that A^^'' < X^^ we obtain 

I - /3ci33 + aCi23p + I - /?C233 + ac223p < A^^^ . (25) 



Next we define dij^ by Cijk = y Xf^ dij^. Then the following matrix is part of a unitary 
matrix if A2^'* 7^ 0: 



Mo 



(^133 C^233 C^333 
C^123 C?223 ^^323 



(2'] 

If A^"^ = 0, we can still choose dijk so that M2 is part of a unitary matrix. Again we perform 
a unitary transformation on the first and second rows as follows: 

^ ^ f a'c?133 + /5'c?i23 a'(i233 + /5'rf223 a'c^333 + /?''i323 

^ \^ -f3'di33 + a'du3 -(3'd233 + "'^^223 -/5'c?333 + a'423 

where 



«' = , ' = , (26) 



AflaP + Afl 



= , ^ = . (27) 

^Af |a|2 + Af 

From the first row of M2 we obtain 

Pcil33 + ;5^cii23r + Rf^233 + ^^^2231' + R433 + ^^23 1 ' < 1 • 

By substituting and ((2Zj), and by recalling the definition of dijk, we find 

|aCi33 + Pcusl"^ + l"C233 + /5c223p + |aC333 + /5c323p < Ag^VP + M^'^I/^T • 

Adding this and and using and pij) we find 

|Cl33p + |Cl23p + |c233p + |c223p + |c333p + |c323p ^ Ag ^ + A3 ^ \a\'^ + Ag • 

Since the left-hand side is equal to x, we obtain x < X2^ + Xf^ by using + = 1. This 
proves inequality (fT7|). □ 



The following result will be useful later. 

Corollary 1. + Ag^'* > A2"^ for all permutations (abc) of (123). 
Proof. Inequality © implies 

2(Af) + Af^) > Af + Xf^ + Af + A^'^) 
> 2A5"^ . 

From this the corollary immediately follows. □ 

V. THE METHOD FOR PROVING INEQUALITIES © © 

All our inequalities have the form 

i^-E(p.A«+rf)+r.Ar)), (28) 

i=l 

where pi, qi and are constants. The function F of the 27 complex variables Cijk with 
constraint ^jjfclQjfcp = 1 is a continuous map from a 53-dimensional unit sphere, which 
is compact, to M. Therefore there must be a minimum value. For all remaining inequalities 
we will show that the minimum of the corresponding function F cannot be negative. 

Consider a three-qutrit state given as in Sec. H] with three diagonal RDMs. Then 
the coefficients Cijk satisfy the orthogonality relations, CijKCiTjK = 0, cijKCjfK = and 
cijkCijk' = 0, where the repeated indices are summed over. The state gives the minimum 
value of F only if the variation of F is nonnegative under any variation 6cijk (not necessarily 
maintaining the orthogonality relations) of the coefficients Cijk- ("A variation" means "a 
first-order variation" throughout this paper.) 

To examine the variation of the eigenvalues X^""^ under a variation of Cijk we use first- 
order perturbation theory in quantum mechanics. An RDM for the ath qutrit, p^"^) , of a 
given quantum state corresponds to the "unperturbed Hamiltonian" , and the variation of 
an RDM under a variation of Cijk to the "perturbation" . The eigenvalues A^^"'* correspond to 
the "energy eigenvalues" . If there are no degenerate eigenvalues, then the variation of A^^"^ 
can be read off from the diagonal elements of 6p[f . We obtain the following result by this 
argument. 

Proposition 3. Suppose that X^^ < X'^^ < Ag""*. Then the variation of the eigenvalues, SXl"'\ 
under a variation Scijk of the coefficients Cijk satisfying jk^^ {^cijkCijk) = are given 
by 

2 E Re {ScukcUk) , 

J,K 

2 E Re {ScijKCijK) , 

I,K 

2 E Re {dcukcHk) ■ 
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Proof. It is enough to note that the right-hand sides give the variation of the diagonal 
elements of the RDMs. □ 

For the function F given by ()28p we define the level L{i,j, k) of the triple [ijk] by 

L{i, j, k) = Pi + qj + Tk . (29) 

Lemma 1. Suppose that the function F of cjjk defined by has its minimum value at 

a state with nondegenerate eigenvalues X^""^ for all a. If Cijk and Ci'j'k' are nonzero for this 
state, then L{i,j, k) = L{i',j', k'). 

Proof. Assume first that i 7^ i', j 7^ j' and k 7^ k' . Consider the variation of coefficients 
Cijk given by 

Scijk = aciijiy , Sci'j'k' = —acijk 
and 6cijK = if [UK] 7^ [ijk], [i'j'k']. Then by Proposition IHl we have the following 
variation of A^"'': 

5AW = 5Af = 5Af = A, 
SX^^ = 6X^^ = 6X^ = -A, 

where A = 2Re {ctCiij'k'CiJk), with all other 6X^^^ vanishing. Hence the variation of F is 

SF = {pi + qj + Tfc - Pi, - qj> - rfc/) A 
= [L{i,j,k)-L{i',j\k')]A. 

It can readily be verified that this formula is valid even if i = z', j = j' or k = k'. Since 
Cijk 7^ and Q/j/fc' 7^ 0, we can make A positive or negative by adjusting the phase of a, 
thus making 6F negative if L{i,j,k) 7^ L{i',j',k'). Therefore, if the function F takes its 
minimum value at a state with nondegenerate eigenvalues, and if Cijk and Ci'j'k' are both 
nonzero, then we must have L{i,j, k) = L{i',j', k'). □ 

Let us define 

(30) 
(31) 
(32) 
(33) 

Our task is to show that the function P^ cannot have a negative minimum for any m. It is 
useful to note that if the quantum state reduces to a two-qutrit system, then Pm > as we 
have mentioned before. We state this fact as a lemma. 

Lemma 2. If Xf^ = 1 and A^"^ = A^"^ = for some a, then Pm > for all m. 

Let L{i,j, k) be the level of [ijk] defined by for the function P^- It can readily be seen 
that 

P^ = Y^L{i,j,k)\ci,k\^ . (34) 

If Pm takes its minimum value at a state with nondegenerate eigenvalues, the level L{i,j, k) 
must be the same for all [ijk] with Cijk 7^ by Lemma [T] This fact can be used to establish 
that Pm cannot have a negative minimum at a state with nondegenerate eigenvalues. 



Pa 


= 2Af 


+ x? 


+ 2Xf 


+ x? 


- 2A« 






= 2Af 


+ Af) 


+ 2Af ) 


+ Ar 


- 2Af ) 


-A(^) 




= 2Af 


+ Af) 


+ 2A?) 


+ xf^ 


- 2A?) 




P7 


= 2X? 


+Ar^ 


+ 2Af 


+ Ar 


- 2A« 
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Proposition 4. // Pm defined by i3(J\} - ^^} for some m has a negative minimum, then the 
state where this minimum occurs must have at least one pair of degenerate eigenvalues, i.e. 
^'f^ — ■^i+i /^'^ some a and i. 
Proof. All Pm have the following general form: 

P„ = 2A;f + + 2Ag' + Ag' - 2A<;> - Ai" . (35) 

Let Z3 be the number in the set {1,2,3} which differs from ii and i2, and similarly for ji 
and ki. There are only four triples [ijk] with negative levels L{i,j, k), which are 

[^iji^i] with L = —2 , 

[kjiki] , [iij2ki] , [iijik2\ with L = -1 . 

Suppose that all eigenvalues are nondegenerate. Then for the state where the negative 
minimum occurs, the level L of the triples [ijk] such that Cijk 7^ must all be —2 or —1 by 

Lemma ^ We must have A3"'' > > for all a. This cannot be the case if there is only 
one nonzero coefficient Cijk- Hence we cannot have L = —2. Suppose L = —1. Then, the 

coefficients Cijk other than Ci^j^ki, Ciij2ki and Ci^j^k2 would be zero. Since X^^^ and must 
be nonzero, we must have is = js = ^3 = 1. This is not satisfied by any Pm- Hence if Pm 
has a negative minimum, then the state where the minimum occurs must have at least one 
pair of degenerate eigenvalues. □ 

Next we establish some general results for degenerate cases. Note that if there are degen- 
erate eigenvalues, we need to make sure that the variation of the off-diagonal elements in 
the degenerate sector vanishes in order to use the perturbation theory argument by taking 
the variation of only the diagonal elements of the RDMs into account. 

Remark. We say that the variation of the RDMs is effectively diagonal if it has no off- 
diagonal elements in the degenerate sector. 

Lemma 3. Suppose that A^^'' = A-^\ 7^ A-,^'' with no other degeneracy at a state where the 
function F defined by / f^) has its minimum. Suppose further that {cijk, Ci+ijk) and Ci'j'k' 
are nonzero, i.e. Cijk ^ or Q+i^fc 7^ 0, and Ci'j'k' 7^ 0, for this state. Then 

L{i,j,k)<L{^,f,k')<Lit + l,j,k). 



Proof. Since = Aj-^\, a unitary transformation in the degenerate sector, i.e. between \i) 
and \i + 1) in the first qutrit, leaves the RDMs diagonal. Hence we can let Cj+ijfc = and 
Cijk 7^ using such a unitary transformation. Then consider the variation given in the proof 
of Lemma Then 

SPi^i-^i = Scijk ■ Cj+ijfc + CijkSciJ^ijk = , 
i.e. the variation of the RDM is effectively diagonal. Define A = 2Re {aCi'j/k'Cijk)- We have 
SXf = 5Af = A and ^aJ,^^ = 6Xf = SXf} = -A as before. However, since Af ^ = X\% we 
must have 5Aj^^\ > ^A^^^-* by definition. Hence we have SX^^^ = A, 6Xf^ = if A > and 
5A2i = 0, 5Af ^ = A if A < 0. Hence 

6F = {pi+i + Qj + Tfc - pif - qj> - rfc/)A if A > , 
6F = {pi + qj + Tfc - pii - qy - rk')A if A < . 

By requiring that 6F > for both cases we find the desired inequalities. □ 
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Remark. Analogous statements for the cases with A^- = A^.,:;^ 7^ A^-, and with A^ = 
^^k+i 7^ ■^k''^ '^^^ readily be established. 

Remark. We often use a unitary transformation in the degenerate sector, which keeps the 
RDM unchanged. A unitary transformation on a qutrit is assumed to be in the degenerate 
sector unless otherwise stated in this and the next two sections. 

Lemma 4. Suppose that xf'^ = A'^\ 7^ A.,^"* and Xf^ = Xfli 7^ X^^ with no other 
degeneracy at a state where the function F has its minimum. Suppose further that 
{cijk, Ci+ijk, Ci,j+i,kj Ci+i,j+i,k) ond Cj/j'fc' orc nonzero for this state. Then 

L{t,j,k)<L{z',f,k')<L{z + l,j + l,k). 



Proof. Since A^- = A^-^^ and A^ = A^-^;^, we can use a unitary transformation to have 
Cj+ijfc = Ci,j+i,k = 0. Then consider the variation given in the proof of Lemma Q and define 
A as in that lemma. Then, as in the previous lemma, the variation of the RDMs is effectively 
diagonal, and we have 5X[^^ = SXf = 0, 5AS+\ = 5Xfl^ = A if A > and 5Xf^ = 5Xf = A, 
(5a5^\ = 5AJ5i = if a < 0. Thus, 

SF = (pi+i + g^+i + - p,> - qj> - rk')A if A > , 
5F = {pi + Qj + Tfc - pi> - qj, - rk>)A if A < . 

By requiring that 6F > for both cases we find the desired inequalities. □ 

Lemma 5. For the function F defined by suppose that Af ^ = aS+\ ^ X^^ and xf = 

^j'+i 7^ ^''^d that there is no other degeneracy. If {cijk,Ci^ijk) o,nd {ci'j'k',Ci'jr+i^k') o,re 
nonzero at a state where the function F is minimized, then 

L{i,j,k) < L{i',j' + l,k') and L{i' , f , k') < L{i + 1, j, k) . 



Proof. We can use unitary transformations to have Cjj+i_fc = Ci'j'^i^k' = and c^fc 7^ 0, 
Ci'j'k' 7^ without changing the RDMs. Then by using the variation given in the proof of 
Lemma n we find 

6F = {pi+i + Qj + rfc - pi> - qj> - rfc/)A if A > , 
6F = {pi + qj + Tfc - pii - qy+i - rk')A if A < . 

By requiring 5F > for both cases we find the desired inequalities. □ 

Lemma 6. Suppose that Af ^ = X^^^ ^ X[}\ Xf = Xf^^ ^ Xf and Xf = A^^ ^ A^^^ at a 
state where the function F defined by i28^} has its minimum. Suppose further that at least 
one of the eight coefficients cuk with I = i,i + 1, J = j, j + 1 and K = k, k + 1 is nonzero 
and that Ciijik' is also nonzero for this state. Then 

L{t, J, k)<L{i',f,k')<L{z + l,j + l,k + l). 
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Proof. Since Xf^ = X\^i, Xf'' = Xf^i and A^^"* = A^^^^, a unitary transformation can be 
used to maximize |cjjfcp. This makes Ci+ijk, Cij+i^k and Cij^k+i vanish. (This argument 
was used to construct a generahzed Schmidt decomposition in Refs. 0, 0]-) Then consider 
the variation given in the proof of Lemma ^ The variation of the RDMs is effectively 
diagonal, and we have SX\^^ = 5Xf = 5Xf^ = 0, 5AS+\ = SXf}^ = SXfl^ = A if A > and 
5Af ) = Xf = 5Af = A, 5X^^, = 5Xfl, = (5Agi = if A < 0. Thus, 

6F = (pi+i + qj+i + Tfc+i - Pi' - qy - ry)^ if A > , 
5F = {pi + Qj + Tfc - pi> - qj> - rfc')A if A < . 

By requiring that 6F > for both cases we find the desired inequalities. □ 
The following proposition will also be useful. 

Proposition 5. Suppose that Xf''^ = Aj-^\ with no other degeneracy at a state where the 
function F defined by has its minimum value. Define the reduced level R^^\j, k) by 

RW(j,k) = qj+rk. (36) 

IfR^'\j,k)^Rmif,k'), then 

CijkCij'k' + Ci+i,jrfcQ+i,j'fc' = , (37) 
i.e. the vectors (cjjfc? Q+ijA;) O'^d (cjj'fc'? Cj+ij'/c') ore orthogonal to each other, for this state. 
Proof. We consider the following variation: 

Scijk = acijiki , 5cj+ijfc = acj+ij/fc' , 
^Cijik' = —acijk , 5ci+ijik' = —aci+ijk 

with all other 6cjjk vanishing. Then we find that the variation of the RDM for the first 
qutrit vanishes. We also find that if we define 

A = Re [a{cijkCij'k' + Q+ijfcQ+ij'fc')] ^ 

then 6Xf = -5Xf = A and ^A^^^ = -SX^ = A. Hence 

6F = (g, +rk- - r,OA = [R'^'\j, k) - R^'\j', k')]A . 

If Eq. (jnZj) does not hold and if R^^\j, k) ^ R^^\j', k'), then the parameter a can be adjusted 
so that 6F < 0. Hence, at a state where F has its minimum value we must have either (jHTj) 

or i?«(j,A;) = i?W(j',A;')- □ 

Remark. The statements obtained from Lemmas 0] and El and Proposition El by qutrit 
permutations are also valid. 

It is useful to treat some cases with only one pair of degenerate eigenvalues in general. 

Lemma 7. Suppose that the function Pm given by ^'^5\l in the proof of Proposition ^ has 
a negative minimum at a state with Xf^ = X\^^ and without any other degeneracy. Then 
Ag' = Ag' = 0. 
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Proof. The only triples with negative levels are [iijiki], [i2jiki], [iij2ki] and [iijik2]. If 
(cjiiifci, Ciii2fci5 Ciijifca) Cij3fc are nonzero, then the preceding lemmas imply that there must 
be some i' (which is not necessarily equal to i) such that L{i',j2, k) is less than L{ii,ji, ki), 
L{ii,j2,ki) or L{ii, ji, k2). This is impossible because L{i',j3,k) > for all i' and k. 
Hence, if Cij^k 7^ for some i and k, then Cj^j^fci = Cjij2fci = Ciijifc2 = 0- Then a unitary 
transformation can be used to have Ci^j^k^ = without changing the RDMs. Thus, we can 
make all coefficients Cijk with L{i,j, k) < vanish. Hence from (|34p we conclude that Pm > 
if > 0. Thus, if the function Pm were to have a negative minimum, then we must have 

(2) (3) 

XjJ = 0. In a similar manner it can be shown that A^^ =0. □ 

Corollary 2. Suppose that the function Pm given by in the proof of Proposition^has a 

negative minimum at a state with X^^^ = X^^^ (Xf^ = A^^"*) and without any other degeneracy. 

Proof. The function Pm can be rewritten as 

p _ 1 , a(1) , \(2) x(2) , ,,(3) , (3) 

Then it is clear that the proof of Lemma [7| applies here after performing a qutrit permuta- 
tions. □ 

Lemma 8. If the function Pm given by / t,'/,5j) has a negative minimum at a state with X[^^ = 
Xll^ and with no other degeneracy, then X^^^ = or xfj = Xf^ = 0. 

Proof To have < we must have Ci^j^ki ^ 0, Ci^j^ki ^ 0, Ci^j^k^ ^ or Ci^j^k2 0. Then 
from the preceding lemmas we find that if Ci,_^jk 7^ 0, then L{i3,j, k) < 0. This is true only 

for j = ji, k = ki. Thus, if we assume that A^-^'' 7^ 0, then Cjgj^fci 7^ 0, and orthogonality 
relations give Cjjji/fci — ^iijiki — 0. If the vectors (Cjij2fci5 ^i2j2k\) &iid {ci-^j-^k2i ^i2jik2) "^^^ linearly 
independent, then there cannot be any other nonzero coefficients of the form Ci-^jk or Cjjjfc 
by Proposition El because R^^\j,k) < for [jk] = [jiki], [jik2] and [j2fci], and R^^\j,k) > 
for all other [jk]. This implies that A^^ = A^^ = 0. If these vectors are linearly dependent, 
a unitary transformation can be used to have Cijj^fc2 = Cjij2fci = 0. Then all coefficients Cijk 
with L{i,j, k) < vanish. Hence Pm > 0. □ 

Corollary 3. // the function Pm given by has a negative minimum at a state with 
'^f^ ~ ('^^k^ ~ ^^k"* ) '"^^^^ other degeneracy, then X^^^ = (Xf^ = 0) or X^^^ = 

Af.' = Ai" = A<f = 0;. 

VI. PROOF OF INEQUALITIES (gD AND dH) 

Now we are in a position to prove inequalities (jH) and The following Proposition, 
which was used in Ref. [1], will be useful. 

Proposition 6. If ii < 12, then \ci^jk\^ + iQjjfcP < X^^^ . 
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Proof. If A'J"* = 0, then Ci^jk = for all [jk]. Since we have \ci^jk\^ < , the desired 

inequahty is satisfied. So we assume 7^ 0. By applying the Cauchy-Schwarz inequality 
to the orthogonality relation Ci^jRCi^jK = 0, we obtain 

I |2| |2<|'\(1)_|„ |2Vx(^)_lr 

\^i\jk\ \^i2jk\ — \'~^i\jk\ )\^i2 \'^i2jk\ )■ 



Thus, 
Hence 



I |2 _|_ I |2 ^ _!2_| |2 _|_ I |2 < 

IQijfel ~r 1^12 jk\ — j--]^^ IQij'fcl ~r \Ci2jk\ _ . 



□ 



Remark. We can prove similarly that iQj^fcp + jcjjjfcP < A^-^"* if ji < j2 and that Ic^fciP + 
IcjjfcaP < A^^ if ki< k2. 

We need the following lemma to prove P4 > 0, i.e. inequality 

Lemma 9. The function P4 has its minimum value at a state with C233 = 0. 

Proof. If Ag^"*, Ag^'' or Ag^"* is degenerate with another eigenvalue, then we can use a unitary 
transformation to have C233 = 0. Hence, if C233 7^ at all states where the minimum of P4 
occurs, then the only possible degeneracies are Xf^ = X''2^ and Xf^ = X''2^ . Suppose that C233 
and another coefficient Cijk are nonzero. Then we must have L{i,j', k') < L{2, 3, 3) for some 
[ij'k'] with X^? = A^^-* and Xf} = Xf\ This is impossible because [233] is the only triple 
with L(2, 3, 3) < —2. Hence there must be a state with C233 = where P4 has its minimum 
value. □ 

Proposition 7. The function P4 is nonnegative. 

Proof. We write P4 in the following form 

P4 = A;^) + 2A(^) + X?^ - A^^) + Af - Af . (38) 
By Proposition IHl and the remark following it we obtain 



C133 


r+ 


C333 r 


< 




C213 


p+ 


C223 ^ 


< 


Af 


C231 


i'+ 


C232P 


< 


Af 



By using these inequalities in we find 

P4 > |ci33p + 1 0333!^ + |C213|^ + 10223^ + 1^231^ + |C232|^ 

+A<'> + A«-Af-Af. 
By letting A^^'* = it |cijfcP, and similarly for A3^'', A3^'* and Xf \ we find 

P4 > |ClllP + |Cii2p + |Ci2l|^ + |Ci22p - 2|c233r 

+ |c31ir + |c312r + |C321P + |c322r • (39) 

By Lemma IHl we can choose a state with C233 = to achieve the minimum of P4. Then, 
Eq. ()39|) shows that this minimum is nonnegative. □ 
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Our next task is to prove P5 > 0. 
Lemma 10. //C133 = C132 = 0, then P5 > 0. 
Proof. The function P5 can be written as 

P5 = X? + 2A« + Xf - + Xf - . (40) 
By Proposition IHl and the remark following it we have 

|c233r+|c333r < Xi^\ 

|cn3p + |ci23r < . 
By using these inequalities in (pn|) we find 

P5 > |C233|'+ |C333|' + |Cll3|'+ |Cl23|' 
+a(^) + A« - Af + Af ) - Af . 

By using the expressions of X'^\ X^^\ X^^^ and Ag^"* as sums of 1%^^ we find 

A > Af)-|ci3l|'-|ci32p-2|Ci33p 

+ |C211|^ + |C212|^ + \C221\^ + \C222\^ 
+ |C311|'+ |C312|' + |C321|'+ |C322|'. 

Since A^^"* > |ci3ip, we have P5 > if C133 = C132 = 0. □ 

Lemma 11. The function P5 is nonnegative if X^^ = X'^\ X^i^ = X^^ or X^2^ = X^r^\ 
(2) (2) 

Proof. If A2 = A3 , then the inequality P5 > is equivalent to 

2Af + Af ) + 2Af ) + Af - 2A« - A^^^ > Q. (41) 

Since Ag^"* > 1/3 > A^^^ we have Xf^ + X^i^ > X^^^ Hence inequality (j4H) will follow from 

A?^ + Af ) + Af ) + Af - X? - A« > . 

This follows from inequality (PJ) by letting (abc) = (123) and noting that A3 > A2 • If 
Af) = , then 

P5 = 2Af + Af ) + 2Af + Af ) - 2X? - A^^ 
= P4 + A« - A« . 

This is nonnegative because P4 > 0. Finally, let A2 = A3 . From Corollary^ we have 

Af + A« > Af . 
By adding this to inequality © with (abc) = (321) we obtain 

2Af + Af ) + 2A« + A« > 2Af + Af = 2Af + Af . 
This is equivalent to the inequality P5 > 0. □ 
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Thus we only need to deal with the cases with A^^"* = X^2 \ ■^2^'' = ■^3^'' ■^f^ = '^2^''- 

Lemma 12. Suppose that the minimum of is negative. Then it does not occur in any of 
the following three cases: 



1. 


A« = A(^) 


with 


no 


other 


degeneracy; 


2. 


A^^) = A^^) 


with 


no 


other 


degeneracy; 


3. 


Af ) = Af 


with 


no 


other 


degeneracy. 



Proof. Case 1. By Lemma |H1 we have either A3 = or A2 = A^^ = 0. The first case is 
impossible. The second case leads to A3^'' = 1. Then, by Lemma El we have P5 > 0. 

Case 2. By Lemma [U we must have A2^'* = A^^"* = 0. This implies that A3^'* = 1. By 
Lemma El we have P5 > 0. 

Case 3. By Corollary |21 we must have Ag^'' = A^^'' = 0. This is impossible. □ 

The following lemma completes the proof of inequality (jSj). 

Lemma 13. Suppose that the minimum of P^ is negative. Then it does not occur in any of 
the following three cases: 



1. 


A« 


= A« 


— H ' 








2. 


A« 


= A« 


and a\ 


= Af 


with no 


other degeneracy; 


3. 


A« 


= A« 


and a\ 


= Af 


with no 


other degeneracy. 



Proof. We will show that there is a state with C133 = C132 = where the function P5 has its 
minimum if it occurs in any of the three cases. Then by Lemma El we conclude that the 
minimum must be nonnegative. 

Case 1. Any unitary transformation on the first qutrit leaves the RDMs diagonal. Con- 
sider the vectors (0133,0233,0333) and (0132,0232,0332). By a unitary transformation on the 
first qutrit we can set 0133 = 0233 = 0132 = 0. 

Case 2. Suppose that (0132,0232) = (0,0) or (0133,0233) = (0,0). Then both 0132 and O133 
can be made zero by a unitary transformation on the first qutrit. Hence, all we need to 
show is that if (0132,0232) 7^ (0,0), then (0133,0233) = at a state where the function P5 is 
minimized. Suppose (0132, O232) 7^ (0, 0) and 033^ 7^ for some k. Then by Lemma (21 we 
have L(3, 3, k) < L{2, 3, 2) = 0, which implies that = 0, i.e. k = 3. Next, suppose that 
(o3ifc,032fc) 7^ (0,0). Then, by LemmaElwe have L(3, l,fc) < L(2,3,2), i.e. < —1. This is 
impossible. Hence (031^,032^) = (0,0) for all k. Thus, 0333 is the only nonzero coefficient of 
the form 03^^. Then by orthogonality relations we have (0133, 0233) = (0, 0). 

Case 3. Note that the triple [133] is not linked to the degenerate sector of the RDMs in the 
sense that aS^^ ^ A^^\ A^^ ^ A^^ and A^^ ^ A^^^ Since L(l,3,3) = -2 and L{t,j,k) > -2 
for all other triples [ijk], if O133 7^ 0, then Cijk = for any other [ijk]. This implies that 
Ag^^ = 0, which is a contradiction. Thus, we can assume that O133 = 0. Suppose that O132 7^ 0. 
If (o33fc,023fc) 7^ (0,0) as well, then by LemmaElwe have L(2,3, A;) < L(l,3,2) = —1, hence 
rfc = 0, i.e. k = 3. Next, if (032^,031^,022^,021^) 7^ (0,0,0,0), then by Lemma H we have 
L(2, 1, A;) < L(l,3,2). This is impossible for any k. Hence 032^ = c^ik = C22k = 02ife = 0. 
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Thus, the only nonzero coefficient of the form c^jk {c2jk) is C333 (C233). Then, by orthogonahty 
we must have C233 = since Ag^^ 7^ 0. Then A2^'* = and hence X^^^ 7^ Ag^'', contradicting the 
assumption. □ 



VII. PROOF OF INEQUALITIES © AND dTj) 

For the case where there are two pairs of degenerate eigenvalues, it is useful to consider 
the variation of the eigenvalues for some variation of the coefficients cjjk in a general setting. 

Proposition 8. Suppose that xf^ = X\^i and Xf^ = X^^^i, and that there is no other 
degeneracy. Suppose further that 

^ijk \ {cijik',Ci+i,j'k') are nonzero. 

Ci+ljk Cj+lj+l,fc / 

Let 

Cj,i+i,fc — Cj+ij,fc — by using unitary transformations, and let 

SCijk = OliCijik' + Q;2Q+1j'A:', = PlCij'k' + f32Ci+l,j'k' , 

SCij/k' = —OTiCijk — , ^Cj+ij/fc/ = —a^Cijk — f32Ci+l,j+l,k , 

where j' 7^ j and j' 7^ J + 1; but k' may or may not be equal to k. Define Oj+i = 
2Re (aj+iCi+ijik'Cijk) and bj^i = 2Re (/3j+/Cj+/j/fc'Ci+ij+i,fc), / = 0, 1. Then it is possible 
to arrange the variation of the coefficients Qj+i^fc and Ci+ijk so that the variation of the 
RDMs is effectively diagonal without making ai, 02, 61 and 62 all vanish. If k ^ k' , the 
variation of the eigenvalues of the RDMs is given as follows: 

5Xf^ = -\a2 -bi\, 5A-^\ = 1^2 - &i| , 

SXf^ = min (fli + a2, bi + 62), ^^f+i = ^lax (ai + 02, 61 + 62) , 
SXy = -ai - 02 - 61 - &2 , SXl = ai + 02 + 61 + 62 , 
6Xf} = -ai - a2 - 61 - 62 • 

If k = k' , then 6X^^^ = SX^j^} = 0, and the variation of the other eigenvalues remains the 
same. 

Proof. By requiring the variation of the RDMs to be effectively diagonal, we find 



Cijk^Ci+lJk + ^Cjj'+l,fcCj+l,j+l,fc — ~^ [Cjj'fc'Q+lj'A:'] ; 



j+i^k + Sci+ijkCi+ij+i,k = . 

Thus, if |cjjfcP 7^ |ci+ij+i,fcP, then we can solve for 6ci^ijk and ^Qj+i^fc to make the vari- 
ation of the RDMs effectively diagonal. If |cjjfcp = |cj+ij+i^fcp, we ffist use phase trans- 
formations to have Cijk = Cj+ij+i^A;. Then any unitary transformation on the ffist qutrit 
can be compensated by a unitary transformation on the second qutrit while maintaining 
the condition Ci+ijk = Cij+i,A,. = 0. Thus we can make Ci+ij/k' vanish as well. Then by 
letting 6ci+ij'k' = Sci+ijk = ^Cj,i+i,A; = 0, the variation of the RDMs can be kept effectively 
diagonal. We have 02 = ^2 = ^ and 02 = &2 = as a consequence but can have ai 7^ 
or 61 7^ by the assumption on the coefficients. (Since we have made Ci+ij'k' vanish, the 
equation 02 = /32 = does not give rise to any further constraint on the variation of the 
eigenvalues.) The calculation of the variation of the eigenvalues is straightforward. □ 
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Remark. This proposition remains valid after a qutrit permutation. 

Now we prove Pq > 0, i.e. inequality ©• Our proof is rather lengthy. 
Lemma 14. //A^^ = X^^\ A^^ = A^^^ or A^^ = A^^^ then Pq > 0. 
Proof. Inequality Q with (abc) = (123) reads 

Af + Af ) + + Af > A^) + A« . (42) 
From inequality we have 

Af + Af > i (Af + aS^) + Af + Af )) 

> ^(a« + aw 



By adding these two inequalities we obtain 

2Xf^ + Af ) + 2A?) + Af > A« + A^^) + A« . 

This is equivalent to inequality (jHI) if A^"*^^ = X^\ Next we recall that A3^^ + X^2^ > X^^ 
(Corollary Q). By adding this to ()42j) we have 

Xf> + Xf + X? + 2Ar) + Xf > 2A^) + A« . 

This is equivalent to inequality (0) if Ag^"* = Ag^"*. Next, if we add the inequality Aa^"* + Ag^"* > 
A2^'* to inequality (jUJ, we have 

2Af + X? + Af + 2Af > 2A« + A« . 

This is equivalent to inequality (0) if Ag^"* = Ag^"*. □ 

The following lemmas will be useful in proving Pq > 0. 

Lemma 15. The function Pq is nonnegative if C231 = C233 = C331 = 0. 

Proof. By Proposition IHl we have |c2iip + |c22iP < Ag^"*. Hence 

Pe = X? - A« + 2A?) + Af) + Af - Af ^ 

> A^ — A2 "* + aP + A2 + X^^ — A^^'* + |c22lP + |c21lP • 

Then, if we substitute A,-^'' = J2jk kijii'P and the similar formulae for A^^"* and A^^'*, the only 
coefficients that contribute negatively to Pg are C231, C233, Cggi. Hence, if these coefficients 
vanish, then Pe > 0. □ 

Lemma 16. The function Pq is nonnegative if C232 = C233 = C332 = 0. 
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Proof. Note first that 

Pe > P6 = 2Af + Af) + 2Af ) + Af - 2A« - A« . 

It can be shown that Pe > under the assumption of this lemma by the argument in the 
proof of the previous lemma with A^^'* ^ A2'^''. □ 

Lemma 17. If C233 = C333 = 0, then Pq>0. 

Proof. Since A?^ + A^^^ > A^^^ by Corollary OQ we have 

Pe > Af + Af ) + 2Af - A(^) - A« 

> \P + Af ) + Af + Af ) - A« - A(^) . 

Using the formulae for the eigenvalues in terms of |cjjfcp in the last expression, we find that 
only |c233p and IC333P have negative coefficients. Hence Pe > if C233 = C333 = 0. □ 

First we treat the cases with only one pair of degenerate eigenvalues. 

Lemma 18. // the function Pg has a negative minimum, then it does not occur at a state 
with any of the following properties: 



1. 


X? 


= A^" 


with 


no 


other 


degeneracy; 


2. 


Af> 


= Af 


with 


no 


other 


degeneracy; 


3. 


Af 


= Af 


with 


no 


other 


degeneracy. 



Proof. Case 1. By Lemma |H1 we have A^ = or Ag = ^2 =0. The latter case leads 
to Ag^'' = Xf^ = Ag^'' = 1, and the inequality is obviously satisfied. Hence we can assume 
that A^^^ = 0, i.e. Cijk = for all [jk]. Then A^^^ = A^^^ = 1/2. If (0231,0331) = (0,0), 
then by letting C233 = by a unitary transformation on the first qutrit and using Lemma IT^ 
we have Pq > 0. If (0231,0331) 7^ (0,0), then all other vectors of the form {c^jk.c^jk) must 
be orthogonal to (0231,0331) by Proposition El Thus, if we let 0331 = by using a unitary 
transformation on the first qutrit, then 02jfc = for all [jk] except [31]. Hence I0231P = 
A^^^ = 1/2, and \ f > I0231P = 1/2, which is a contradiction. 

Case 2. Corollary |2 implies that A^^^ = A2^'* = 0, and by Lemma |21 we have Pg > 0. 

Case 3. Assume that O233 7^ 0. Then, since there is no [ij] satisfying L(i,j, 3) = — 1 
except [ij] = [23] and since the triple [233] is not connected to the degenerate sector, the 
O233 is the only nonzero coefficient of the form Cij^. Then, by orthogonality relations we 
have 0231 = 0232 = 0. If (oiji,Oij2) 7^ (0,0), by LemmalHwe must have L{i,j, 1) < L(2,3,3), 
i.e. Pi + Qj > —1. Hence [ij] = [23], which has been rejected, [21] or [33]. Then A2'^^ = 0, 
and we have Pg > by Lemma |21 Thus, we can assume that O233 = 0. If (0231,0232) = 
(0,0) [(0331,0332) = (0,0)], then we can have O331 = [0231 = 0] in addition by a unitary 
transformation on the third qutrit. Then Pg > by Lemma ITHl Assume that (0231, O232) and 
(0331,0332) are both nonzero. These vectors are orthogonal to each other by Proposition 
applied to the case with degeneracy in the third qutrit. Then the only other possible nonzero 
vector of the form (04^1,0,^2) is (0211,0212) by the same proposition. For Ojjs to be nonzero 
we must have L{i,j, 3) < min [1/(2, 3, 2), L(2, 1, 2), L(3, 3, 2)] = 1. All nonzero coefficients of 
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the form q^s must have the same value of L{i,j, 3). In order not to have X^^ = (and hence 
Pe ^ by Lemma 121) we must have L{i,j,3) = 1 and the possibly nonzero coefficients of 
the form Cjjs are C223, C313 and C133. (We mean, strictly speaking, that Cijs is nonzero only 
if [ijS] is in the set {[223], [313], [133]}.) Then C233 = C333 = 0, and by Lemma IT7I we have 
Pe > 0. □ 

Remark. We abbreviate the phrase "possibly nonzero coefficient" as PNC below. 

The following lemma is used in the most degenerate case. 
Lemma 19. //C231 = C211 = C233 = C213 = 0, then Pe > 0. 
Proof. First we note inequahty (jHI) with (abc) = (123): 

A^^) + X? + Xf + Af - A?' - A^^) > . 
The inequality Pe > follows from this if 



Af + Af - A?) > , 



By using the expression for X'^^ in terms of \cijk\'^ we find that this inequality is satisfied if 
C211 = C213 = C231 = C233 = 0, as required. □ 

Now we treat the most degenerate case. 

Proposition 9. // the minimum value of Pg is negative, then it does not occur at a state 
with X'^^ = Xf\ X^i^ = Ag^"*, Af"* = X^2^ and with no other degeneracy. 

Proof. Note that we can make C233, C231 and C211 vanish by successively using unitary trans- 
formations on the first, third and second qutrits. Then the only PNC with a negative level 
is C331. We will show that if C331 7^ 0, then C213 = 0. Then by Lemma ITUl we have Pq > 0. 
Thus, we assume that 

f C331 C332 \ / C323 C313 

\ C231 C232 / \ C223 C213 

are nonzero and derive a contradiction. 

Suppose that the first matrix has rank less than two. Then after making C233 vanish 
by using a unitary transformation on the first qutrit, we can make C231 and C331 vanish 
simultaneously by a unitary transformation on the third qutrit. Then, by Lemma El we 
have Pe > 0. Thus, we can assume that the first matrix has rank two. We set C231 = C332 = 
by unitary transformations on the first and third qutrits and let C313 = by using a unitary 
transformation on the second qutrit. Then, the above matrices become 

° ] and ° 

C232 / V '^223 C213 

with C331 7^ and C232 7^ 0. Let us consider the following variation: 

5C323 = -«lC33i - a;2C232 , 

^C213 = — /5iC33i — /32C232 , 

5C223 = -71C33I - 72C232 , 

Scssi = aTC323 + /5lC213 + 7rC223 , 

^C232 = Ct2C323 + /32C213 + 72C223- 
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By requiring that the variation of the RDMs be effectively diagonal we have 



C331'5C231 + 5C332C^ + (5c3i3C^ = -5 [C323C^] , 

^C3i3C323 = —5 [C213C223] ; 
C232'5C231 + 5c332Ciir = . 

These equations can be consistently solved for 5c23i, ^0332 and 6C313 provided that C323 7^ 
and |c33ip 7^ |c232p- (Notice that the variations 5c23i, ^0332 and ^0313 do not contribute to 
the variation of the diagonal elements.) If C323 = 0, then we can use a unitary transformation 
on the second qutrit to eliminate C223 as well. Then 5[c323C223] = ^[£2130223] = if we set 
7i = 72 = 0. Then we can simply let ^0332 = 5c3i3 = 5c23i = 0. If IC331P = |c232p, then 
we can use a phase transformation to have C331 = C232. Then any unitary transformation in 
the first qutrit can be compensated by a unitary transformation in the third qutrit to keep 
C331 and C232 unchanged and keep the condition C231 = C332 = 0. This means that we can 
diagonalize the second matrix as well, making C223 vanish in addition. Then we again have 
^[c323C223] = 5[c2i3ci^] = if we let 7i = 72 = 0. Then we can let ^0332 = 6C313 = 6C231 = 0. 

Let us introduce the following definitions: = 2Re (aiC(*^C323), 6, = 2Re (/?iC^*^C2i3) and 
Cj = 2Re (7iC'^*^C223), where C^^^ = C331 and C*-^-' = C232. The variation of the function 
can be given as follows: 

6Pe = 5AS^) - + 6Xf - 6Xf + 6Xf - 6Xf . 

Assume first that C323 7^ and |c33ip 7^ |c232p and that the variation of the RDMs is made 
effectively diagonal. Then, we have 

SPq = |a2 - 61 - ci I - min (ai + a2 + Ci + C2, 61 + 62) 
-(ai + 61 + ci + 02 + 62 + C2) 
+ \a2 + 62 + C2 - ai - 61 - ci| . 

If Oi can be made nonzero, then so can 02, and vice versa, and similarly for bi, 62, Ci and C2 
because C331 and C232 are both nonzero. Hence we can always have 02 + 62 + C2 = ai + &i + Ci 
by adjusting a^. Pi and 7j, i = 1,2. We can also choose all coefficients to be nonnegative. 
Then, 

SPq < |a2 - &i - Ci I - (ai + 61 + Ci + 02 + &2 + C2) . 

If ai and 02 can be made nonzero, then by choosing the others to be zero (by letting 
Pi = 'ji = 0, i = 1, 2) we find 6Pq < —ai < 0. If ai = 02 = 0, then 6Pq < —62 — C2. Since 
62 or C2 can be made positive — otherwise C323 = C223 = C213 = — we have 6Pq < 0. If 
C323 = 0, then we have ai = 02 = and Ci = C2 = to make the variation of the RDMs 
effectively diagonal, but we can still make 61 and 62 positive. Then we can have 6Pq < 0. 
Similarly, if IC331P = |c232p, we must have Ci = C2 = 0, but we can still have either ai 7^ 0, 
02 7^ or 61 7^ 0, 62 7^ (or both). Then we can have 6Pq < 0. 

Thus, if C331 7^ and C213 7^ at the same time, we can always lower the value of Pq 
if it is negative. Hence, if C331 and C213 are both nonzero, the function Pg does not have a 
negative minimum. □ 

Finally we deal with the cases with two pairs of degenerate eigenvalues to complete the 
proof of inequality (jHl). 
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Lemma 20. Suppose that the minimum of the function Pq is negative. Then it does not 
occur at a state with X^^ = X^^^ and Xf'^ = X^^ and without any other degeneracy. 

Proof. If (c23i,C33i) = (0,0), then we can use a unitary transformation on the first qutrit 
to make C233 vanish. Then Pq > hj Lemma IT31 Thus, we can assume that (0231,0331) 7^ 
(0,0). If (0233,0333) = 0, then we have Pq > hj Lemma IT71 Thus, we can also assume 
(c233,0333) 7^ 0. By Proposition El we conclude that the vectors (0231,0331) and (0233,0333) are 
orthogonal to each other. By the same proposition (0232, O332) must be orthogonal to both 
of these vectors. This implies that (0232, O332) = (0, 0). Then by a unitary transformation on 
the first qutrit we can have O233 = 0, and we conclude that Pe ^ by Lemma ITO □ 

Lemma 21. Suppose that the minimum of the function Pg is negative. Then it does not 
occur at a state with X^^ = X^^^ and xf^ = X^^ and without any other degeneracy. 

Proof. Suppose that (0231,0331) = (0,0). Then we can have O233 = while maintaining 
C231 = O331 = by a unitary transformation on the first qutrit. Then Pq > hj Lemma IT31 
Next, suppose O213 = O313 = O223 = O323 = 0. First we use a unitary transformation on the 
first qutrit to let O233 = 0. Then a unitary transformation on the third qutrit can be used 
to have 0231 = 0. Since 0213 = 0223 = 0, by Proposition IHl we have 

1^ |2 I 1^ |2 < x(3) _ x(3) I ,2 I |2 
IC331I + IC333I ^ ^3 — ^3 ~ |C213| ~ |C223| • 

We use this together with |o2iip + I0221P < A2^^ to have 

Pq = A« - A« + - A?) + 2Ar + Xf^ 
> aW - A« - Xf> + 2Af 

+ |c21iP + |022l|^ + I0331P + 1 0333!^ + |02l3p + 10223^ • 

By substituting the expression of Aj"'' in terms of loyfcp we find 

P6>-2(|o23ir + |o233r)=0. 

Thus, it is enough to show that if (0231, 0331, 0232, C332) 7^ (0,0,0,0), then (o223,C323) = 
(C213,C313) = (0,0). 
If the matrix 

/ O231 0331 
\ O232 0332 

has rank less than two, then we can set O231 = O331 = by a unitary transformation on the 
third qutrit. Then a unitary transformation on the first qutrit can be used to let O233 = 0. 
Then Pq > 0. Therefore we can assume that the matrix is regular. 

First we use Proposition |H1 after a qutrit permutation (123) (132) with [ijk] = [213] 
and [i'j'k'] = [132] to show that (0223,0323) = (0,0). By considering the variation as given 
in Proposition IHl we find 

SPq = 5AW - 5A« + 5Af - 5Af + 5Af - 5Xf^ 

= |a2 - ^1 1 - max (ai + 02, &i + 62) - 3 min (ai + 02, 61 + 62) • 

This can be made negative unless ai + a2 = or 61 + 62 = 0. Since the matrix in is 
regular, we can have Oj 7^ if 6, 7^ and vice versa for both i = 1 and 2. Thus, unless 
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Oi = 02 = &1 = &2 = 0, we can have both oi + 02 and bi + 62 nonzero. Hence SPq can 
be made negative. Thus, unless (0223,0323) = (0,0), Pq cannot have a minimum. Next we 
change k' from 2 to 1 to show that (0213,0313) = (0,0). The only difference here is that we 
have 6X'^^ = instead of X^2^ = —ai — 02 — 61 — &2- Thus we find 

SPe = 1^2 - &i| - 2min (ai + 02, h + 62) • 

Again, this can be made negative unless 01 + 02 = or 61 + 62 = 0. Hence by the 
same argument as above we conclude that the function Pg cannot have a minimum un- 
less (0213,0313) = (0,0). □ 

Lemma 22. Suppose that the minimum of the function Pg is negative. Then it does not 
occur at a state with Xf'^ = X^2^ and Xf^ = X^^ and without any other degeneracy. 

Proof. Assume that 0233 7^ 0. Then by lemmas in Sec. |3 we find that the PNCs other 
than 0233 are (0231,0232), (0211,0212,0221,0222) and (0331,0332). By a unitary transformation 
we can let O212 = O221 = 0. By orthogonality relations 0233023^ = 0, i = 1,2, we have 
0231 = 0232 = 0. Hence the PNCs are 0233,0211,0222,0331 and 0332. If 0332 = 0, then I0331P > 
I0222P + |o233p + I0211P but |o233p > |o33iP and I0222P > |033iP- These imply that 0331 = 
C233 = O211 = O222 = 0, which is a contradiction. Hence O332 7^ and by an orthogonality 
relation we have O331 = 0. Then 

Pe = "~|C233P ~ ^211^ + 2|0222p + I0332P > 

because I0332P > I0233P + I0211P + |o222p- 

Next suppose that 0233 = 0. Suppose further that (0231,0232) = (0,0). Then, we can use 
a unitary transformation on the third qutrit to have O331 = 0. Then Pg > by Lemma 1131 
Hence we assume that (0231, O232) 7^ (0, 0). Let us use a unitary transformation to let O232 = 
for convenience. Note that L(2, 3, 1) = —2 and L(2, 3, 2) = 0. Since L(l, 3, 3) = 1, we must 
have 0133 = by LemmaEl For (0113, 0123) to be nonzero, we must have L(l, 1, 3) < L(2, 3, 2), 
which is false, by Lemma El Hence (0113,0123) = (0,0). Similarly, for (0313,0323) 7^ (0,0) we 
must have L(3, 1,3) < L(2,3,2), which is false. Hence (0313,0323) = (0,0). Thus, the PNCs 
of the form Ojj3 are 0213, 0223 and 0333. (We have already shown that 0233 = 0.) 

By applying Proposition El to the case with degeneracy in the third qutrit, we find that 
(0131,0132) and (0331,0332) are mutually orthogonal and are both orthogonal to (0231,0232). 
Since we have let O232 = 0, we obtain 0131 = O331 = and either 0332 or 0132 must vanish. 
However, if 0332 = 0, then we have Pg > by Lemma ^1 Hence we can assume that O332 7^ 
and 0i32 = 0. Next, we use Proposition |H1 after a qutrit permutation (123) (321) with 
[ijk] = [111] and k' = 2. Then we find 

6Pq = 3(ai + 02 + 61 + ^2) + 2|a2 - bi \ + max (ai + 02, 61 + 62) • 

By making all and 6, nonpositive, we obtain 

5P(i < 3ai + 02 + &i + 362 • 

The right-hand side can be made negative unless Oi = 02 = 61 = 62 = 0. Hence Om = O112 = 
0121 = 0122 = 0. 
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Next we apply Proposition |S1 with [ijk] = [113]. The only difference is that we have 
6X^1^ = instead of 6X^i^ = ai + a2 + 6i + &2- Hence 

6Pe = 2(ai + a2 + h + 62) + 2|a2 - &i| + max (ai + 02, 61 + 62) ■ 

If Oi or 62 is nonzero, then 6Pq can be made negative by making the nonzero one large 
and negative. If 02 and bi are both nonzero, then we can let 02 = &i < ai = 62 = 
to make 6Pq < 0. If 6Pq > and if C231 7^ and C232 = 0, then the PNCs among c^jk, 
j, k = 1,2, is C322 after using a unitary transformation in the second qutrit. If we undo this 
transformation, we have C312 7^ as well. 

Next let [ijk] = [112] in Proposition |S1 Then ^A^^"* = SX'^^ = and all other variations of 
the eigenvalues are the same. Hence 

5Pe = ai + 02 + 61 + &2 + 2|a2 - &i| + max (ai + 02, h + 62) • 

If tti or 62 is nonzero, then by letting the nonzero one be negative and letting a2 = fei = 
we have 

SPe<ai + b2<0. 

If ai = 62 = 0, then 

SPe = 02 + fei + 2|a2 -bi\+ max (02, 61) . 
If both a2 and bi are nonzero, then by letting 02 = &i and letting them be negative, we find 

5P6 < 02 + &i < . 

By the same argument as in the k = 3 case, we have only C212 and C222 among C2jk, j, k = 1,2, 
possibly nonzero if 6Pq > 0. Here we can use a unitary transformation on the second qutrit 
to let C212 = 0. 

Thus, the following coefficients are the PNCs: 

C332, C322, C312, C222, C213, C223; C333, C231. 

Recall that C332 7^ 0. Hence, by orthogonality relations we have C322 = C312 = 0. We may 
assume that C213 7^ because otherwise A^^"* = A2^'* = and -Pe ^ by Lemma El Then by 
the orthogonality relation C213C223 = 0, we have C223 = 0. From C332C333 = we have C333 = 0. 
Then C233 = C333 = 0, and by Lemma IT7I we have Pe > 0. □ 

Now we turn our attention to inequality ((7j), i.e. Pj > 0. 

Lemma 23. The function Pj is nonnegative if X^i^ = X^2 \ ^^i^ — ■^2'^ '^'^ -^2^^ — -^s^"*- 
Proof. If A^^^ = A2^\ then P7 > is equivalent to 

2Xf'> + A?) + 2Af + Xf > 2X? + X^i^ . 

Since A3^'' > 1/3 > Af'* and A^^'' > 0, it is enough to show that 

Af ) + Af + Xf + Xf > aS^) + A« . 
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This follows from inequality Q since 

Ar > ^1 • If ^1 = ^2 > then P7 > is equivalent to 

2Af + Af ) + 2Af + Af > 2A(^) + A« , 

which follows from inequality Finally, if X'^^ = Xl^\ then Xf'^ = 1 — 2A2^''. Hence 

P, = 2- 3A?) + 2Af + Af - 2A« - A^^^ 
= 2 - 3A(^) + Xf - Af ) - A« + aS^) . 

Since Ag^"* and X^^ are both less than or equal to 1/2, we have P7 > 0. □ 
Thus, the only degeneracies we need to consider are A2^^ = Ag^^ Xf^ = X'^^ and X^2^ = Xf\ 
Lemma 24. If C231 = C221 = C331 = 0, then Pj > 0. 
Proof. Writing Pj as 

Pr = aS^) - A« + X? - Xf + 2Xf + Af 

and noting that |c232p + IC233P < Xl^\ which follows from Proposition IHl and the remark 
following it, we have 

P7 = ^^2^ + '^3^'' + -^1 ~" -^2 "* + -^1 "* ~ -^3 "* + |c232p + |c233p • 

Then, by expressing the eigenvalues in terms of Icy^p we find that the negative terms are 
-2|c23iP - |c22iP - |c33iP- Heucc, if C231 = C221 = C331 = 0, then P7 > 0. □ 

Lemma 25. If X^^^ = A^^^ and Xf^ = A^^^ = A^ , then P7 > 0. 

Proof. Consider the vectors (C231, C232, C233) and (C331, C332, C333). We can use a unitary trans- 
formation on the third qutrit to make C231 and C331 both vanish. Then we use a unitary 
transformation on the first qutrit to let C221 = 0. Note that the condition (C231, C331) = (0, 0) 
is preserved. Hence, C231 = C331 = C221 = and P7 > by the previous lemma. □ 

Lemma 26. Suppose that the minimum of the function P7 is negative. Then it does not 
occur at a state satisfying any of the following conditions: 



1. 


A« 


= A^" 


with 


no 


other 


degeneracy; 


2. 


Af> 


= Af 


with 


no 


other 


degeneracy; 


3. 


Af 


= Af 


with 


no 


other 


degeneracy. 



Proof. Case 1. If (C231, C331) = (0, 0), then a unitary transformation on the first qutrit can be 
used to have C221 = 0. Then we have P7 > by Lemma YIM Next let (C231, C331) 7^ (0, 0). By 
LemmalHlwe have A^^^ = or A^^'' = X^^^ = 0. The latter is impossible. Hence A^^^ = 0. Then 
X'^^ = Ag^^ = 1/2. Use a unitary transformation to let C331 = 0. Then, since all {c2jk,C3jk) 
must be orthogonal to (0231,0331) for [jk] 7^ [31] by Proposition we have C2jk = except 
for C231. This implies that |c23ip = Ag^"* = 1/2. Hence A^'^^ > 1/2. This is impossible. 
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Case 2. Assume that (0231,0232) 7^ (0,0). If Qjs 7^ for some [ij], then we must have 
L{i,j,S) < L(2,3,2) = —1. This is impossible because = 2. If = for all [ij], 

then Xf^ = 0, which is impossible. Hence C231 = C232 = 0. If Cij^ 7^ and if (0221,0222) 
or (0331,0332) is not (0,0), then we must have L{i,j,3) < L(2,2,2) = L(3,3,2) = 0. Since 
r3 = 2, we must have pi = -2, qj = 0, i.e. [ij] = [23]. Thus, if (0221,0222) or (0331,0332) is 
not (0, 0), then the only nonzero coefficient of the form 0,^-3 is O233. If the vectors (0221, O222) 
and (0331,0332) are linearly dependent, then a unitary transformation can be used to have 
0221 = 0331 = 0. Then by Lemma ED we have P7 > 0. If they are linearly independent, by 
Proposition El applied to the case with degeneracy in the third qutrit we conclude that all 
the other vectors (ojji,Ojj2) must vanish since [ij] = [22] and [33] are the only pairs with 
R^^\i,j) = pi + Qj = -1. Then the PNCs are 

C221; C222; C331 , O332, O233. 

We can let O221 = by a unitary transformation on the third qutrit. Then by an orthogonality 
relation we find O332 = 0. [Note that the vectors (0221,0222) and (0331,0332) are hnearly 
independent.] Hence, O222, C233 and O331 are the PNCs. Then I0331 

1 2 = = A^^^ = 1/3 and 

■^2^^ = |c222p + |c233p = 2/3. This is impossible. Hence, we must have O231 = 0221 = O331 = 
and P7 > by Lemma 

Case 3. Assume that O231 7^ 0. If {cij2,Cijs) 7^ (0,0), then L(z,j, 2) < L(2,3, 1) = —2. 

This is impossible. Since Cijs = for all [ij] would imply that A3 = 0, we conclude 
that 0231 = 0. If 0221 or 0331 is nonzero and if (ojj2,Ojj3) 7^ (0,0), then we must have 
L{i,j,2) < L(2,2,l) = L(3,3,l) = —1. This is satisfied only by (0232,0233). Then, by an 
orthogonality relation we have O232 = or 0233 = 0. Hence, either A2 =0, which yields 
P7 > by Lemma or A3 = 0, which is impossible. Thus, 0231 = 0221 = 0331 = and 
P7 > by Lemma El □ 

Lemma 27. Suppose that the minimum of the function P7 is negative. Then it does not 
occur at a state with A^ ■* = A2^'' = A3^'* and without any other degeneracy. 

Proof. First we establish the following result, which is similar to Proposition El if 

{ciji, Cij2, Cijs) and (oj/j/i, Oj'j/2, Oi/j/3) are nonzero, then either R^^\i,j) = R^^\i'j') or these 
vectors are orthogonal to each other. This can be shown as follows. Consider the variation 

SCijK = CiCi'j'K , SCi'j'K = —CiCijK foT all K 

and ScijK = for all other coefficients. We note first that the variation of the third RDM 
vanishes. The variation of P7 is found to be 

3 

k=l 

Hence, if 6Pj > for all a, we must have either R^^\i,j) = R^'^\i' ,j') or Y^k=i (^i'j'kCiJk = 0. 
We note that R^^\2,3) = -2, R^^\2,2) = R^^\3,3) = -1, and R^^\iJ) > for all other 

If 023fc = for all k, then we can have 0221 = 0331 = by a unitary transformation on 
the third qutrit and have P7 > by Lemma EH Assume that (0233,0232,0231) is nonzero. 
Let O232 = O231 = by a unitary transformation on the third qutrit. Then, if 033^ = for 
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all k, we can use a third-qutrit unitary transformation between |1) and |2) to have C221 = 
while keeping C232 = C231 = C331 = 0. Then we have P7 > by Lemma 1^ So we assume 
that C33A; 7^ for some k. By orthogonality of (C231, C232, C233) and (C331, C332, C333) and by the 
condition C231 = C232 = we have C333 = 0. We choose C332 7^ and C331 = 0. If C221 = 
for the vector (C221, C222, C223), then P7 > by Lemma So let C221 7^ 0. Then, all the 
other vectors (cjji, Cjj2, Cij^) are orthogonal to three linearly independent vectors. Hence they 
must all vanish. Then, the PNCs are C233, C332 and C221 since C222 = due to the relation 
C222Ciir = 0. Then, IC332P > |c233p + |c22iP and |c233p > IC332P imply that IC332P = |c233p, 

^3) 

C221 = 0. Hence X\ =0, which contradicts the degeneracy assumption. □ 
The following lemma completes the proof of inequality ((Zj). 

Lemma 28. Suppose that the minimum of Pj is negative. Then it does not occur in either 
of the following two cases: 

1. \^2^ = X^^^ and X^2^ = X^^'' with no other degeneracy; 

2. X^2^ = X^^^ and Xf^ = X^2^ with no other degeneracy. 

Proof. Case 1. If (0231,0331) = (0,0), then we can use a unitary transformation on the first 
qutrit to have C221 = 0, and use Lemma EH to conclude that P7 > 0. Suppose (0231,0331) 7^ 
(0,0). If (cij2,cij3) 7^ (0,0), then we must have L(l,j, 2) < L(3,3, 1) = —1, which is 
impossible. Hence (cij2, C1J3) = (0,0) for all j. Next let (02^2, £3^2, C2j3, £3^3) 7^ (0,0,0,0). 
Then by considering the variation given in Proposition |S1 after a qutrit permutation (123) — *• 
(132) with [ijk] ^ [22j] and j' = 1, we find 

5P, = 5AS^) - 5X^^ + 6X? - 5Af + 5Af - 6Xf^ 

= \a2-bi\+ max (ai + 02, h + 62) + (4 - j){ai + a2 + h + 62) • (44) 

If j = 1 or 2, then by choosing a, and 6, to be nonpositive, we have 

SPj <\a2-bi\+ 2(ai + 02 + &i + 62) < 

if some of a, and 6, is nonzero. For j = 3, if Oj and 6, are nonpositive, then 

5Pf = \a2-bi\+ 2 max (ai + 02, h + 62) + min (ai + 02, hi + 62) • 

If ai 7^ or 62 7^ or if both 02 and bi are nonzero, then we can have 6P7 < 0. If 02 is the 
only nonzero coefficient, then we find 6 Pi = |a2| + 02 = 0, and similarly for bi. The cases 
02 7^ and 61 7^ are related by unitary transformations. By choosing the latter, we find 

that C233 = C332 = C333 = 0. Since we have already concluded that there are no coefficients 

(3) 

of the form Cij^ other than C233 and C333, we have A3 =0, which is a contradiction. 

Case 2. It is enough to show that (0231,0331) = (0,0) for the same reason as in case 1. 
Assume that 

(2:2:)^(oo)- («) 

If this matrix has rank one, then we can use a unitary transformation on the third qutrit 
to have C231 = C331 = 0. Then, we can have C221 = using a unitary transformation 
on the first qutrit and conclude that -P7 > by Lemma |^ Therefore we may assume 
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that the matrix ()45|) has rank two. Suppose further that cijs ^ 0. Then we must have 
L(l, j, 3) < L(3,3,2) = 0, which is impossible. Hence cya = for all j. Next, suppose that 
(c2j3,C3j3) 7^ (0,0). Then by letting (abc) = (132) and [ijk] [213] and j' = 3, k' ^ j in 
Proposition IHl we find 

SPj = \a2 - bi \ - min (ai + 02, h + 62) - (4 - j)(ai + 02 + 61 + ^2) • 

This formula will be identical with ()44|) if we let —ai and bi ~bi. Hence 5Pj can be 

made negative if (0213,0313) or (0223,0323) is nonzero. For j = 3, then we can have SP^ < 

unless 02 or bi is the only possibly nonzero one among Oj and 6j. Since the matrix (j45p has 

rank two, we must have both 0231 and O332 nonzero after diagonalizing it. Hence, if Oj can 

be made nonzero, then so can 6j, and vice versa. Hence we can conclude that unless a, and 

bi, i = 1,2, are all zero, we can have 5P7 < 0. Hence we must have Cijs = for all [ij], i.e. 
(3) 

X's'' = 0, which is a contradiction. □ 



VIII. CONSTRUCTION OF SOME STATES 

We have shown that the inequalities listed in Theorem^ are necessary for a state to exist 
with a given E-point. Our next task is to show that these inequalities guarantee that there 
is a quantum state having these eigenvalues. The following elementary fact will be useful. 

Lemma 29. Suppose that the angles 6 and ip are constrained by 

u sin 29 — V sin 2(/9 = , 

where u and v are real constants, and that there are no other constraints on them. Then, 
9 + if can take any value and the range of the function 

f = u cos 29 + V cos 2{p 

is given by \\u\ — |f|| < |/| < \u\ + |f |. 

Proof. The constraint equation can be written as 

(m — v) sm{9 + if) cos(6' — (p) + {u + v) cos(6' + if) sm.{9 — </)) = . 

For any value oi 9 + ip, we can find a value of 6* — so that this equation is satisfied. Hence, 
the combination 9 + is unconstrained. Note that 

= ('UCOS26' + t'cos2(y9)^ + ('usin26' — wsin2(y9)^ 
= + v'^ + 2uv cos 2(6^ + v?) • 

This shows that |/| ranges between (m + f )^ and {u — v)'^. Hence, the range of / is given 
by — |f II < I/I < \u\ + \v\. (Since / —f for 9 ^ 9 + 7i/2 and ip + 7i/2, the range is 
invariant under f —>■ —f.) □ 

First we construct a convex subset of S. 
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Proposition 10. There is a state with any E-point in the convex set Sq bounded by hyper- 
planes with the following comer points: 

[B,B,B], [B,A,A], [A,B,B], [0,B,B], [A, A, A], 

[A5,|M]' [A HI, OH] 

and those obtained by qutrit permutations from them. 

Proof. First consider the following PNCs: 0333,0322,0232,0223,0311,0131,0113. Suppose that we 
add C122, C212 and C221 to this list. Then the only off-diagonal elements in the RDMs that 
become nonzero are the ones connecting |1) and |3) in each qutrit. Thus, each RDM can be 
diagonalized by a simple unitary transformation. Motivated by this observation, we consider 
the state with the PNCs given as follows. First let 

f &333 \ ^ f acos9i\ f 6331 \ ^ ( -9 sin Oi \ 
V ^133 ) y a sin 6*1 y ' \^ 6131 J \ gcosdi J ' 

( &313 -hsinO, \ ( 6311 \^(fcofieA 

\bn3 J \ hcosOi J ' \bni J {fsmO,)' 

where a, h, g and / are real constants. Next we define 

f disk \ ^ f cosifi - sin \ f h^k \ 
\ dilk J V / V ^^1*; / ' 

where [ik] = [11], [13], [31] or [33]. Finally we let 

( ^'^^ ^ ^ cos xi - sin Xi\ f dij3 \ 
\Ciji J \smxi cosxi J \diji J 

with [ij] = [11], [13], [31] or [33]. We also let 

O322 = b cos 62 , O122 = b sin 62 , O232 = cos (p2 , O212 = sin ip2 , 
0223 = d cos X2 , 0221 = d sin X2 , 

where 6, and d are real constants. Then we find a'^ + b'^ + + d'^ + f'^ + g'^ + = 1 and 





= c' 


+ d' 


A? 


= b' 


+ d' 


xr 


= b' 


+ c' 



and 

A^^) - aS'^ = b^ cos 202 + (a' + h^) cos 26^ , 

Af ) - Xf^ = cos 2ip2 + (a^ + 9^ - f^ - h'^) cos 2^i , 

Xf - Xf = d^ cos 2x2 + (a' + /i' - - g^) cos 2xi . 

The orthogonahty relations are 

6^ sin 2^2 + (a^ + - - h^) sin 2^i = , 
sin 2(/?2 + (a^ + 9^ - f^ - h^) sin 2v7i = , 
sin 2x2 + (a^ + - f^ - g'^) sin 2xi = . 
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We require that 



2 I 2 

a +g 



g-h^ > 0, 
f-h' > 0, 



Then, by Lemma 1221 we conclude that Ag"^ — A^"'' have the ranges given by the following 
inequahties: 



I 2 I <-2 2 

\a +f -g 

I 2 I 2 p2 

\a +g -f 



^2 I 



< A 



(1) 

3 

.(2) 



A« < a' 

.(2) 



f-g'-h' + b' 



Af ) <a^ + h^-f-g^ + d\ 

(a) , (a) 



(46) 
(47) 
(48) 



The equations for Ag and the inequalities for Ag — A]^ can be satisfied for all corner points 
listed here. We give a^, 6^ etc. that are nonzero for each corner point: 



[B,B,B] 
[B,A,A] 
[A,B,B] 
[0,B,B] 
[A, A, A] 
[0,A,A] 
[0,0,0] 

[AB, III] 

L ' "4 4' 4 4 2j 



1/4, 
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c' = dP = 1/6, f 
6^ = 1/3, = d^ = 1/6, 
1/3, 6^ = 1/12, ^ = d? = 1/4, f 



= 1/12, 
= 1/12, 



1/2, 



b^ 
b^ 
1, 

(f = 1/3, 
1/2, = 



1/3, 
= 1/4, 



f = 1/6, 
1/4. 



The E-points obtained by qutrit permutations can be realized because this construction is 
invariant under qutrit permutations. If two sets of eigenvalues A^^"^'^ and A-"^'^'^ are given by 



6^ = bj, etc. and 



aXl""^' + (1 



fj, 6^ = bjj, etc., respectively, then any set of eigenvalues 
a)X'f^^^ with < a < 1 is given by = aaj + (1 — a)a//. 



of the form A^^'^'* 

b^ = ab'j + (1 — a)b'jj, etc. with inequalities (pH|) - (jlH|l being satisfied. Hence the set of states 
constructed here forms a convex set in the space of E-points. □ 



It will be useful later to construct states whose E-points are on the boundary of inequality 
((T}. This inequality with (abc) = (123) reads 

Pi = Af ) + Xf + Xf + Af - aS^) - a(^) > . 



for all [ijk] except 

3. We also find that unless A2^^ = Ag^"*, we must have C332 = C331 = 0. 



From the proof of Theorem |21 we find that to have Pi = we need Cijk 
those with z = 3 or j 
Thus, a generic state satisfying the equality Pi = has only the following PNCs: 

Cl31, C231, Ci32, C232, C113, C213, C123, C223, C333. (49) 

Notice that for all these coefficients Cijk we have L{i,j, k) = 0. The orthogonality relations 
are as follows: 

0, (50) 
0, (51) 
. (52) 



C131C231 + C132C232 + C113C213 + C123C223 
C113C123 + C213C223 
C131C132 + C231C232 
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The last two constraints are satisfied by letting 

cii3 = acos6', C2i3 = asm9, Ci23 = —bsm9, C223 = bcos9, 
Cm = d cos ip, C231 = -d sin ip, C132 = / sin ip, C232 = / cos ip, 

where a, b, d, f, theta and (p are real numbers. Let us write the eigenvalues A^-"'' as A^-"'' for 
later purposes. Then 

Af^ = a', K^^^ = b\ Kf = d\ Kf = f 

and 

A« + aS^) = Af + kf^ + Af + Ap) , (53) 
A^') - aS'^ = (A^'^ - Af ^) cos 26 + (Af - Af ^) cos 2^ . (54) 

The constraint (jSUj) reads 

(A^'^ - kf^) sin 25 = (Af - kf) sin2v3 . 

Hence, by Lemma 1^ the range of A« - A« is given by 

|(Af - Af )) - (A^^) - Af ))| < A^^) - A« < Af - A^^ + kf - kf . (55) 

In summary, there is a state with any E-point satisfying (jSHj) and and A2'^^ > A^"'* 

for all a. Notice that it is not necessary for the eigenvalues to satisfy Ag*^"* > Ag""*, or even 
Ag"'' > A^"''. It is clear that the states constructed here form a convex set in the space of 
I3-points. It can readily be seen that the corner points of this boundary hyperplane listed 
in the proof of Proposition [T] all satisfy and, of course, Hence all E-points on 

the boundary hyperplane (j^Hj) have corresponding states. This construction of states can be 
slightly generalized to prove the following lemma, which will be useful later. 

Lemma 30. There is a state with any E-point in the simplex with the following corner 
points: 

[0,0,0] [B,0,B] [B,B,0] [A,0,A] [A,A,0] [i?,0||,Oi|], [B,A,A]. 



Proof. Let us add C322 to the list of PNCs ()49|) . Then we find that no new orthogonality 

relation will be introduced. Let us define A^^"* to be the sum over J and K of |cjji^p, where 

CijK appear in the list ()49|1 . and similarly for A^-^"* and A^^'*. Then the eigenvalues A-"^ can 
be expressed as follows: 

AP = A[ ^ A2 = A2 \ A3 = A3 -* + |c322p , 
X?=k?, X?^=k?^ + \C322\\ 

X? = kf\ Ar=Ar + |c322p. (56) 

The preceding construction shows that if A^"'' are given in this manner and if A^"'' satisfy 
and (jSSl), then there is a state with the eigenvalues X\'^\ Furthermore, if we require 
that A2"^ > A^"'' for all a, then the set of E-points satisfying all these conditions is convex. 
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All points except [B, A, A] satisfy the conditions given here with C322 = since they are 
on the boundary (j^Hj) with A^-'^^ = A^'^^ for all i and a. Thus, all we need to do is find A^^'^'* 
[which satisfy (j^Hj) and and C322 such that the eigenvalues [B, A, A] are expressed as 
in We can do so with |c322p = 1/6, A? = A? = 1/3, A^^^ = 1/6, Af^ = Af^ = 0, 

□ 



a(2) = A^''^ = 1/3 (and with A^"^ = A^^^ = 1/2 as a resuh). 



,(3) 



,{3) 



We will need the states given by the following PNCs: 

C132 = cos 9 , C232 = —a sin 9, 

C112 = r sin a sin 9 , C212 = r sin a cos 9, 

C122 = r cos a sin 6* , C222 = r cos a cos 9, 

Cm = fsm9, C231 = / cos 9, 

C113 = p sin a cos 9 — q cos a cos ip , 

C213 = —p sin a sin 6^ — g cos a sin , 

C123 = P cos a cos 6* + g sin a cos ip , 

C223 = — P cos a sin 6* + g sin a simp , C333 = (7, 

where a, /, (7, p, q, r, 9, (p and a are real numbers. (Note that the PNCs will be those for 
the boundary Pi = if r = although the parametrization is slightly different.) The only 
nontrivial orthogonality relations are 



p^)sm29 = —g^ sin 2^9. 



(r^ + p^ — q^) sin 2a 



2pq cos{9 + ip) cos 2a . 



(57) 
(58) 



The eigenvalues are given by 



(1 



A 



(1 



(2 



A 



(2 



(3 



Ai' 



A 



(3 



(a^ + p^) cos^ 6* + (r^ + /^) sin^ 6* + g^ cos^ , 
(a^ + p^) sin^ 6^ + (r^ + f) cos^ 6^ + g^ sin^ , 

/ 2 I 2\ • 2 I 2 2 

(r + p )sm a + q cos a 



(^2 



pq sin 2a cos(6' + (p) 
p^) cos^ a + g^ sin^ a + pq sin 2a cos(6' + y?) 



2 I 2 



Let us introduce the variable = + q^. Then we have 

A« + A« = a^ + r^ + f + s^ 
Af + Af = 



2 I 2 
r + s , 



Af 



2 I 2 



Al' 



(3) 



Notice that these equations can be solved for a?, r^, p and 
be satisfied by letting 



(59) 
(60) 
(61) 
(62) 

Then note that Eq. (|58|) can 



sin 2a 



cos 2a 



2pq cos{9 + (p) 



(r^ + p2 _ _|_ 4p2g2 cos2(6' + ip) 

^2 _|_ p2 _ g,2 

-^Z (r^ + — g^)^ + 4p2g2 cos2(6' + ip) 
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unless + — = pq cos{6 + (/?) = 0. The eigenvalues A^^^ and A2^'' can be given as 



.(2) 



(2) 



A 



(2) 



2 
1 

2 L 



r 



+ + - a/ (r2 + p2 _ g2)2 _^ 4p2g2 cos2(6' + if) 



r 



+ + q^ + ^(r2 + p2 _ ^2)2 + 4p2^2 cos2(^ + y,) 



These equations are valid even if + — = cos (6' + v^) =0. 
Suppose that a^, r^, and are fixed at some values satisfying 



s2_^2 > g 



f-2r^ > 0. 



We first show that A^^^ can take any value between and (r^ + s^)/2. By the first part of 
Lemma l29l we find that the combination 6 + !f can vary freely. Hence A^^"* can take any value 
as long as 



1 



+ S^- ^(r2 + p2 _ ^2)2 + 4p2^2j < ^(2) < ^ [^2 ^ ^2 _ |^2 ^ ^2 _ ^2|j 



By varying q^ from to (r^ + s^)/2, which is possible by the assumption that < s^, we 
can see that A^^^ varies from to (r^ + s^)/2. 

However, let us impose the condition < A^^'' or, equivalently, X'^^ < s^. We will show, 



under this condition, that Ag"*^^ — A^^"* can take any value between — p — 2r^ — X)^' + A^' 
and - p - 2r^ + A^^ - Af^|. Note first that 



(2)^ X(2)| 



Af^Af = g2[r2 + p2gin2(^^^)] 



(63) 



This equation to have a solution for 6' + if and only if 

< Xfhf'^ < q\r^ + p") = q\Xf^ + A^ - q^) . 



These inequalties are satisfied if A^^"* < q^ < Ag^''. Assuming that A^^'' > r'^, we have 
g2 > > and p2 = ^2 _ ^2 > ^2 _ ^^2) ^ ^(2) _ ^2 ^ q_ rj.^^^ ^^^^ ^ ^.^2^^^ ^ 



(2) 



,{2) 



function of with A^^"* fixed as 



sin^(6' + v?) 



Af ^Af - gV^ 

p2g2 



By substituting this in 

|A?)-A«r 



[(r^ + /2 - - cos 2^ - g^ cos 2(/?] ^ 
+ [{r^ + f -a^ - V^) sin 26 + g^ sin 2^?] ' 

/ 2 I ^2 2 2 2\2 

(r +/ -a -g) 

+4(r2 + /2 _ - p2)g2 gj^2^^ ^ ^ 



we obtain 



|A«-A«P 



2\2 



(r^ + /2 - - p2 - g 

\(2)\(2) _ „2„2 
I ^ / 2 I <-2 2 2\ ^1 ^2 y ' 

+4(r +/ -a ) 5 
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If we substitute q'^ = X[ and = A2 , then the right-hand side becomes \a'^ — p — 2r^ + 
■^2^^ — ^ p and ~ P — 2r^ — + Af ^ respectively. Hence, by continuity, | Ag^^ — 
takes any value between these values as ranges from to A^^^. If A^ ^ = r^, then Eq. ^ 
can be solved by letting q^ = = Ag^"*. Then we have = 0, and sin^(^ + (p) can take any 
value in [0, 1]. Hence, X^2^ — Xf'^ can take any value between — p — 2r^ — X^2^ + Xf'^ \ and 
|a2-/2-2r2 + A^^^ + Af^|. This range contains the desired range. The result obtained here 
can be summarized as follows. 

Lemma 31. Suppose that 







+ X?^ 


+ Af) 


+ xr 


-A«-A«>0, 


a' 


- A« 


+ A« 


~X? 


-Af 


-Af)>0, 




- A« 


+ A« 


-xf 


-Af 


> 0. 



Let p = Xfl Ifr^ < Af) < A^ < s\ - p - 2r^ > and 

\a' -p- 2r' - A?) + Af ^ | < A« - A« < a' - p - 2r' + Xf^ - Xf^ , 

then there is a state with these eigenvalues. The set of E-points satisfying these conditions 
is convex. 

Proof. The first three conditions guarantee that Eqs. (jKHjl . (|Hnjl and can be solved for 
a^, and s^. Convexity is obvious from the form of the conditions imposed. □ 

The following fact, emphasized by Bravyi [3], is useful in constructing states. 

Lemma 32. Consider the set of states such that Cijk 7^ only if the triple [ijk] is in a set 
E. Suppose that if [ijk] and [i'j'k'] are in E and if [ijk] 7^ [i'j'k'], then at most one of 
the equations i = i' , j = j' and k = k' holds. (This implies that there are no nontrivial 
orthogonality relations.) Then the set of E-points for these states is convex. 

Proof. Let A^"''^ and A^"''^^ be the eigenvalues of the states with the specified nonzero coeffi- 
cients. Let Aijk = |cjjjtp and write Aij^ = Al-^^ and Al^j. for these states. Then, for any A["'* 
given by 

A}'^) = aXf^' + (1 - «)A;'^)'' 
with < a < 1 for all a and /, a corresponding state is obtained by letting 

A^,k = «4fc + (1 - a)All^ 
for all [ijk]. □ 

IX. CONSTRUCTION OF THE SET 5 

We have constructed states with some E-points in the convex set S in Proposition IIUI and 
Lemma IHUl We will start constructing states with their E-points in outlying subregions of 
S. First we divide the region S into two: one region Si that satisfies 

P3 ^ 2Xf + Xf + 2Af ) + Xf - 2A(^) - A^^) < (64) 
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and the other which satisfies Ps > 0- The only corner point of S with Pg < is [A, iii, Qif ] 
which has Pg = —1/4. The corner points satisfying Pg = are 

[0,0,0], [B,B,0], [A,A,0], [A,0,A], 

[A, B, A] , [Oif , B, 0||] , [A, III] , [A, B, |||] . (65) 
Using this observation, we show that the set 5*1 consists of three simphces. 

Lemma 33. The subset of S satisfying Pg < 0, i.e the convex set Si, consists of the three 
simplices Ci, C2 and C3 with the following five E-points being corner points of all three 
simplices: 

[0,0,0], [B,B,0], [A,A,0], [A,B,\\l], [A,\\\,Q\1] , 
and with the additional corner points for each simplex given by 



Ci 
C2 
C3 



[A,B,A], [0||,P,0||]; 
[A,0,A], [A,B,A]- 
[A,0,A] , [A, III, III] . 



Proof. Since the E-point [A, ;j|^,0||] is the only E-point satisfying Pg < 0, the boundary 
hyperplanes of Si are Pg = and those of 5" that contain |^|,0^|]. The latter are 

aJ^) = 0, A^^^ = A^^\ Af^ = A^^\ Af^ = 0, © with {aba) = (231), © with {aba) = (132), © 
with (aba) = (132) and ((Tj) with (aba) = (132). It is enough to show that Ci U C2 U C3 is 
bounded by these hyperplanes. 

Let us define the following functions: 

Pg ^ 2AS^) + A^^) + 2Af) + Af - 2Af ) - X?\ 

Pio - A?) - Xr - A?) . (66) 

Then the boundary hyperplanes of the simplex Ci can be given as follows, where the E-point 
listed with each boundary hyperplane is the one which is not on the hyperplane: 

[0,0,0] : ©,(a6c) = (132); [B,B,0] : A^'^ = 0; 
[A,A,0] : Af) = A?); [A5,|M] ■ ^ = 0- 
[A Hi Oil] : ^8 = 0; [A,B,A] : ®, (a6c) = (132); 
[0||,P,0||] : P9 = 0. 



Those for Co are 



Those for are 



[0,0,0] : A« = A?); [B,B,0] : A« = 0; 
[A,A,0] : aS^) = A(^); [A,B,lll] : A^^ = 0; 
[AHI'OH] ■■ Ps = 0; [A,0,A] : P9 = 0; 

[A,B,A] : Pio = 0. (67) 

[0,0,0] : X^^ = X?; [B,B,0] : A^'^ = 0; 

[A, A, O] : Af ) = Af ; [A, P, |||] : @, (abc) = (231); 

[AiiiOif] : ^8 = 0; [A,0,A] : ®, (a6c) = (132); 

Hi' Hi] • = 0. (68) 
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The E-point [0||, B, 0||] is the only corner point of Ci, C2 or C3 with Pg < 0, and all other 
corner points of Ci satisfy Pg = 0. Hence all E-points in Ci satisfy Pg < 0. The E-point 
[A, III, |||] is the only E-point of Ci, C2 or C3 with Pio > 0, and all other corner points of 
C3 satisfy Pio = 0. Hence all E-points in C3 satisfy Pio > 0. All corner points of C2 satisfy 
Pg = = except [A, O, A] with Pg = 1 and Pio = and [A, B, A] with Pg = and 
Pio = —1/3. Hence all E-points in C2 satisfy Pg > and Pio < 0. Thus, the intersections 
Ci n C2 and C2 n C3 are hyperplanes Pg = and Pio = 0, respectively. The corner points of 
these boundary hyperplanes are 

Pg : [A, P, A] , [A, Oif] , [A, B, |||] , [A, A, O] , [P, P, O] , [O, O, O] ; 
Pio : [A, O, A] , [A, lii, Oi|] , [A, 5, Mi] ' ' O] , [O, O, O] . 

Thus, the simplices Ci, C2 and C3 are glued together at these boundaries, and all other 
boundaries are boundaries of the convex set 5*1. Hence S*! = Ci U C2 U C3. □ 

Now we construct states with their E-points in the simplices Ci, C2 and C3. 

Lemma 34. There is a state with any E-point in Si = Ci U C2 U C3. 

Proof. Since C2 C 5*0, there is a state with any E-point in C2 due to Proposition [THl Next 
we construct a state with any E-point in the simplex Ci. From the list ()fi7|l we see that the 
boundary obtained by omitting [A, P, A\ is 

\f + + \f + Af ) - A^^) - A« = 0. (69) 

This suggests that the states may be constructed by using Lemma ED after letting 
{\^i\\^2^) (A^^\A^^^) and This is allowed since the only in- 

equalities among A^-"'' that we need for this lemma to hold are A^"^ > aS"^ for a = 1, 2. For 
an E-point satisfying equality (j69|) . we have, using the notation of this lemma. 





= 


+Ar^ 


+ xf 


+ xf 


-A« 




= 0, 




= A« 


+ A« 




-Af 


-Af) 


= Ar 


>o. 


(2) 
3 


= A?^ 


+ A« 
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and f'^ = X[ . Hence we have < Ag < A3 < s^. We also have 

a2_/2_2^2^_^(3)__^(3)^Q_ 

The possible range of A^^^ - A^^^ IS 

|,(3) ,(3) ,{2) ,{2)| >(1) \W ^ x(3), x(2) ,{2) 

1^2 ~ ^1 ~ ^3 + ^2 I — ^3 ~ ^2 ^ ^2 ~ ^1 + ^3 ~ ^2 • 

These are equivalent to the following three inequalities if Eq. (j69|) is satisfied: inequality (0) 
with (abc) = (231), inequality (jH)) with (abc) = (132) and inequality ((Tj) with (abc) = (132). 
Hence these are all satisfied by the corner points satisfying (I69|) . For [A, P,yl], we find 
r2 = 1/6, a2 = 1/3, s^-Xf^ = 1/6, anda2-/2-2r2 = 0. Thus, |2r2+/2-a2±Af ^-A^^l = 0, 
and the condition 

\a' -f- 2r' - Xf + Af I < A« - A^^^ < a' - f - 2r' + A^ - xf 
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is satisfied by Therefore, by Lemma IHTl all E-points in Ci have corresponding 

states. 

The construction of states with their E-points in C3 is similar. We note that the boundary 
of C3 obtained by omitting [A, B, |||] is 

A« + A« + \f + \f - ) - Af = . (70) 

This suggests that we can apply Lemma |^ with Ag^"* ^ \^^\ Ag^"* ^ Ag^'' and then with 
(1) (2) (3) (1). For an E-point satisfying (f7n|) we have, using the notation of this 
lemma. 
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+ A« 


+ Af) 


+ xf 


-Af) 
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-Af) 
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-A« 
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+ Ar 
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-Ar 


= Af) 



and p = \i . So these are all nonnegative on the E-points satisfying (fTUj). For [A,B, |||] 
we find = 1/6, = 0, = 1/3, — Ag^"* = 0. For the E-points satisfying (j7U|) we have 
a'^-p- 2r2 = A^^^ - aS^^ > 0. For [A, B, |||] we have - P - 2r^ = 0. Then, for C3 to 
be realized by the states constructed here it is sufficient to have 

\2P + p-a' + Af - Af)| < Af - X? < -2r' - P + a' + Af - X?^ (71) 

for all corner points. This is satisfied by [A, B, |||] since A2^'' — Af = Ag^"* — A^ ^ = 0. For 
the E-points satisfying (fTOI) . Equation (f7T|) is equivalent to inequality (jH) with (abc) = (123), 
inequality (jH)) with (abc) = (132) and 

2A« + A« + 2Af ) + Af > 2AS^) + A^^ . 

This inequality can readily be verified by using inequality Q with (abc) = (231) and A3^'' > 
1/3 > Xf\ □ 

Thus, we have constructed a state with any E-point in S with an additional condition 
Pg < 0. Since the intersection of S and the hyperplane Ps = is the 5-dimensional convex 
set with corner points (fHHji . the convex subset of S with the condition Pg > are those of 
S excluding [A, j\^,Ojj]. By symmetry we can construct a state with any E-point in S 
satisfying 

2Xi'^ + Xi'^ + 2Af ^ + A^'^ - 2A^'') - A^^^ < (72) 

for any (abc). Thus, our next task is to construct a state with any E-point in S with 
additional conditions 

2Af ) + X? + 2Xf + Af) - 2A^"^ - A^'^^ > (73) 

for all {abc). This set of E-points is the convex set with all the corner points inherited from 
5* except for 0||] and those obtained by qutrit permutations from it. Let us denote 

this set by 5*2. From the construction of Si it is clear that inequalities ©, ©) (© and ^ 
are redundant once we impose (j73p . We give a more straightforward proof of this fact. 
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Proposition 11. Inequalities (0), (0^ and ^ follow from inequalities (Op, ^ and 

Proof. Inequalities © and (|7j) follow immediately from inequality (fTSjl by noting that A2^'' > 
A^"* and X^/^ > respectively. Inequality (0) is equivalent to 

2a('^) + A^'^) < 2X^^^ + Xf^ + 2Xf + A^'^ . 

By adding this and inequality (f7H|) and dividing by three we obtain inequality ©• Inequality 
(H)) is equivalent to 

2A^'^) + aS'^) < 2Xf + Af ) + 2A^^^ + A^ . 
On the other hand, inequality (fTSj) is equivalent to 

2X^^^ + A^") < 2Af ^ + A^') + 2AS^) + A^ . 
By adding these inequalities together and dividing by three we obtain inequality (j2I)- □ 

We observe that among the corner points of 5*2 (in fact of S) the E-point [5, 0||, 0||] is 
the only E-point which violates the inequality 

P, ^ 2X? + Xf + 2Af ) + Xf - 2X? - XS^ > . 

It has Pg = —1/3. The E^points which satisfy Pg = are 

[0,0,0], [A,A,0], [A,0,A], [B,0,B], 
[B,B,0], [B,A,A], [B,A,lll], [P,|M,A]. 

These E-points are the corner points of the convex subset of 5*2 with the additional condition 
Pg < as we will show below. 

Lemma 35. Let S3 be the convex subset of S2 obtained by imposing an additional inequality 
Pq < 0. Then Ss = D1UD2UD3 where Di, D2 and D3 are simplices. The following E-points 
are the corner points of all three simplices: 

[P,Oi|,0||], [0,0,0], [B,B,0], [B,0,B], [B,A,A]. 

The additional corner points of each simplex are as follows: 

D, : [A,A,0], [B,A,lll]; 
D2 : [A,A,0], [A,0,A]; 

D3 : [A,0,A], [B,lllA]. (74) 

Proof. Since the E-point [P,0||,0||] is the only E-point satisfying Pg < 0, the boundary 
hyperp lanes of S3 are Pg = and the boundary hyperp lanes of S2 that contain [P, ||, 0||]. 

The latter are A^^^ = aJ\ A^^^ = X'i\ Xf^ = 0, Af^ = 0, ^ with {aba) = (231), with 
(abc) = (321), (P) with {aba) = (123), ® with (aba) = (123) and ® with [aba) = (132). It 
is enough to show that Pi U P2 U D3 is bounded by these hyperplanes. We define 

r, - x(l) x(2) , x(2) , .{3) ,(3) 

^ _ x(l) x(l) , x(2) , (2) , (3) , (3) 
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First we note that all E-points in Di satisfy Qi < and Q2 > 0, all E-points in D2 satisfy 
Qi > and Q2 > 0, and all E-points in satisfy Qi > and Q2 < 0. Let us list 
the boundaries of each simplex together with the E-point omitted to obtain each boundary 
5-simplex: 

D, : [i?,Oi|,0||] : P9 = 0; [0,0,0] : A^^^ = A^^^; 

[B,0,B] : dZSD, (a6c) = (321); [B,B,0] : A^^ = 0; 

[B,A,A] : ©,(a6c) = (132); [A,A,0] : aS'^ = Af^; 
[B,A,lll] : Qi = 0; 

D2 : [i?,0||,Oi|] : P9 = 0; [0,0,0] : A^^ = A^^^; 

[B,0,B] : Af) = 0; [5,5,0] : aS'^ = 0; 
[B,A,A] : dH), (a6c) = (123); [A ^, O] : Q2 = 0; 
[AO, A] : gi = 0; 

D3 : [S,0if,0i|] : P9 = 0; [0,0,0] : A^^^ = A^^^; 

[B,0,B] : Af^ = 0; [P,P,0] : (HBD, (a6c) = (231); 
[P,AA : (0),(a&c) = (123); [A,0,A : \^^^ = 
[B,lll,A] :Q2 = 0; 

Thus, the simplices Di and D2 are glued together at the common boundary Qi = with 
corner points [P,0||,Oi|], [0,0,0], [P,0,P], [P,P,0], [B,A,A] and [A, AO], and Da 
and D3 at = with corner points [B, Oif, Oi|], [O, O, O], [5, O, P], [B, B, O], [B, A, A] 
and [A, 0,A]. The other boundaries are all boundary hyperplanes of S3. Hence, 5*3 = 
P>i U ^2 U D3. □ 

Next we construct states which have E-points in this region. 
Lemma 36. There is a state with any E-point in S3. 

Proof. Note that states with their E-points in D2 have been constructed by Lemma EDI 
In order to construct states with their E-points in Di we consider the following PNCs: 
C333, C113, C223, C121, C231, C322, Ci32. Since there are no nontrivial orthogonality relations, the 
set of 15-points obtained from these coefficients is convex by Lemma EH Hence it suffices 
to construct the corner points of Di using this set. This can be done as follows (with 

Aijk = |Cijfc| ) • 

[P,Oi|,0||] : A32 = A23 = A33 = 1/3; 

[0,0,0] : A33 = l; 

[B,0,B] : A33 = A31 = A32 = 1/3; 

[B, B, O] : A33 = ^113 = v4223 = 1/3; 

[B, A, A] : Ai32 = A223 = 1/3, A22 = A33 = 1/6; 

[A AO] : A33 = A23 = l/2; 

[B,A, : A121 = An2 = 1/6, A33 = A23 = 1/3. 
The simplex D3 is obtained from Di by the qutrit permutation (2) ^ (3). □ 
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Thus, we have constructed states with their E-points in the convex subset of 5*2 satisfying 
Pg < 0. By symmetry we can construct states with their E-points in 5*2 if they satisfy 

for some (abc). Let us denote the convex subset of 5*2 satisfying P^""^^^ > for all (abc) by 5*4. 
Our task is now reduced to constructing states with their E-points in 5*4. The corner points 
of S'4 are those of S excluding [A, j\^,0^^] and [i?,0||,0^|] and those obtained by qutrit 
permutations from them. The following lemma, which allows us to reduce the convex set S'4 
further by removing [A, A, A] from the corner points, can be proved by using Proposition 11111 

Lemma 37. The convex subset of S4 satisfying 

Q3 = 1 + X? - A^^) + X? - + Af ) - Af < (75) 
is a simplex with the following corner points: 

[A, A, A], [0,A,A], [A,0,A], [AA,0], [B,A,A], [A,B,A], [A,A,B]. 
There is a state with any E-point in this simplex. 

Proof. It is straightforward to check that > for any corner point of S'4 not listed here. 
We have = —1/2 for [A, A, A] and Qs = for the other E-points in the list. The simplex 
with the corner points given here is bounded by the hyperplane Q3 = and the boundary 
hyperplanes of S4 that contain [A, A, A], which are aS"^ = and A^"^ = \f \ a = 1,2,3. 
Hence, this simplex is the subset of S4 with the condition (fTSj) . All E-points in this simplex 
can be realized by quantum states by Proposition ITUl since it is a subset of Sq. □ 

Thus, our task is further reduced to constructing the convex set S5 whose corner points are 
those of S4 excluding [A, A, A]. The following proposition almost accomplishes this task. 

Proposition 12. There is a state with any E-point in the convex set S^^ bounded by hy- 
perplanes with the same corner points as S5 excluding [0,A,A]. 

Proof. Consider the following set of PNCs: 

{C3335 C113, Ci3i, C223) C2325 C211, C321, C312, C122} • 

Note that this set is invariant under the qutrit permutation (2) ^ (3). Note also that 
there are no nontrivial orthogonality relations to be satisfied. Hence, the set of the E-points 
realized by the states considered here is convex by Lemma 1221 Therefore all we need to do 
is find the values of Aijk = \cijk\'^ for each corner point of S5 except for [O, A, A]. This can 
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be done as follows: 
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Mi' HI] • ^211 = ^232 = ^23 = 1/6, A333 = 1/2; 
A, iif] : A223 = ^122 = ^321 = 1/6, A333 = 1/2; 
[A, B, III] : ^211 = A312 = 1/6, A333 = ^223 = 1/3; 
[III, 5, A] : A333 = A312 = 1/3, A223 = ^122 = 1/6; 
[B, III, A] : A122 = Aii3 = 1/6, A333 = ^32 = 1/3. 

The E-points which are obtained from these by the qutrit permutation (2) ^ (3) can also 
be obtained by letting Ayk — ^ Aikj in this list. □ 

There are two convex sets S^^ and S^^ where all corner points of 6*5 are realized except 
[A, O, A] and except [y4,A, O], respectively, which are constructed similarly. Therefore all 
simplices whose corner points are also corner points of S5 can be realized except those with 
all three E-points [0,^4,^4], [A, 0,^4] and [A, A, O] as corner points. This observation can 
be used to finish the construction of states with their E-points in 6*5. 

Proposition 13. There is a state with any E-point in S5. 

Proof. Any E-point X in the set S5 can be represented as 

X = ^a,n, (76) 

s 

where Yg are the corner points of and where the are nonnegative and satisfy '^s = 1- 
If the coefficients for the corner points [A, |||, |||], [|||,y4, |||] and [|||, |||,y4] are 
zero, then the corresponding E-point is realized by a state by Proposition EB Also, if the 
coefficient ag for [A, yl,0], [A, 0,^4] or [0,yl, A] vanishes, then, the corresponding E-point 
is realized by Proposition and the remark following it. 

Suppose that the coefficient for [A, |||, |||] is a^i^ and that for [0,^4,^4] is a'^'^ for an 
E-point X in Note that 

2 111411] + A, A] = [A, A, O] + [A, O, A] + [O, B, B] . 

Suppose that > 2a^\ Then 

[A^mm+^'10,AA\ = (««-2a(^))[A,|||,|||] 

+4^ {[A, A, O] + [A, O, A] + [O, B, B]} . 
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Thus, the E-point X can be reexpressed in such a way that the coefficient for [O, A, A\ 
vanishes. On the other hand, if < 2a^2 \ then we can write 

"!'^ \\l -M] + [O, A, A] = - laf [O, A, A] 

{[A, A, O] + [A, O, A] + [O, fi, 5]} . 

Hence, the E-point X can be reexpressed in such a way that the coefficient for [A, Hf, Hf] 
vanishes. By symmetry the same conclusions can be made for the case where the coefficients 
for both [|||,y4, |||] and [A, 0,^4] are nonzero or where those for both [|||, |||, A] and 
[A, A, O] are nonzero. 

Let the coefficients for [|||,y4, |||] and [|||)|||)^] be af^ and af\ respectively, 
and those for [A, 0,^4] and [A, A, O] be and , respectively, for the E-point X. If 
> 202"'' for some a, then the E-point X can be reexpressed in such a way that one of 
the coefficients ctg"^ vanishes. Then by Proposition El we can conclude that E-point X has 
a corresponding state. If of* < 2a^^ for all a, then by successively applying the argument 
above to each qutrit, we can make vanish for all a. (Note that the coefficients a'^'' 
never decrease in this process.) Then by Proposition the E-point X has a corresponding 
state. □ 



X. CONCLUDING REMARK 

It is surprising that there is a very simple condition on the RDMs for the n-qubit pure 
quantum state for all n. The analogous condition for the qutrit system has turned out 
to be more complicated as shown in this paper. Nevertheless, the result. Theorem has 
some structure, and it is probable that there is a proof simpler than the one given here. It 
is natural to investigate the corresponding problem for the n-qudit pure quantum system 
for all n. There are two features in our result which may be true in the general case. One is 
that the set of the allowed E-points is convex. The other is that the boundaries of this set 
are hyperplanes in the space of E-points. We conclude this paper by a proposition which 
suggests that the boundaries are indeed hyperplanes in the n-qudit case. 

An n-qudit quantum state is given by 

1^) = Ci,i,...ijii) (S) \i2) (S) ■ ■ ■ ^ \in) (77) 

l<ii,i2,---,in<d 

with Yliii2 i \^hi2-iJ\^ — where |1), |2),...,|(i) form an orthonormal basis. As before we 

assume that the RDMs are diagonal. Let X^f' be the ith. smallest eigenvalue of the RDM for 
the ath qudit. We regard the set of these eigenvalues as a point in the n{d — l)-dimensional 
Euclidean space parametrized by A^-"'' ,l<i<(i — l,l<a<n and call this set an E-point 
as before. It is convenient to define the lattice E-point Li-^i^-.in the E-point given by 
Xll^ = 1 and, result, Aj-"^ = if i ^ ia- Note that these E-points do not correspond to 
a physical state unless h = i2 = ■ ■ ■ = in = d because they do not satisfy the conditions 
A-"^ < A^"^ for i < j. For a given set of ra-tuples we define to be the convex set whose 
elements A are the E-points given as 

A = ^ ] ^iii2---in^iii2---in ) 
[iii2---in]&N 
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where < Aiia-i™ < 1 and E[ni2-i„]e7V 
proving our proposition. 



1. The following lemma is useful in 



Lemma 38. Consider a state of the form \71\) with the eigenvalues Aj-°\ 



1, d, a 



l,...,n, such that Q^j2...j„ 7^ only if [iii2---in\ ^ N. Then the corresponding E-point 
A = (Af ^) IS m Vn- 

Proof. An El-point A = (A-"^) in V^v is given by 



A 



(a) 



Ul,j2,-,jn)eN 



where < Aj^j^,„j^ < 1 and ^^.^ 
the state (f77j) if we make the identification Aj^j 



A- ■ 

i2...in]€N ^1112. --t 



1. These are the formulae for A-"^ for 
= Ic- ■ • P □ 



Proposition 14. Let U he an open set in the d{n—l)- dimensional space of E-points. Suppose 
that the eigenvalues are not degenerate in U. Suppose further that an E-point A = (A'""*) has 
a corresponding state if and only if f{X) > 0, where f is a function with a continuous and 
nonzero gradient and that the hypersurface /(A) = is connected. Then this hypersurface is 
part of the hyperplane Vn for some set of n-tuples N . 

Proof. Suppose that a state with A^^'^'* = A^^'^'* is in U and on the boundary so that /(A) = 0, 

where A = (A^^"^), and that both Cj^j2---jn ^^"^ Q'^i'j-'-ijj are nonzero for this state. We consider 
the variation 



5 a 



l\%2---ln 



acjiji...ji , 

'■l'2 'n ' 



with all other 5ci^i^...i^ vanishing. Then the variation of the function / is given by 

n d—1 



df 
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df df 
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■ 2Re(aCi'^i^...i/^Q,i2...i„) , 



A=A 



where we have defined df /dX^f' = for all a, because there are no degenerate eigenvalues. 
For 6f to be nonnegative for all a we must have 



df 



t(dX 
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df 



A=A 



a=l 
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(78) 



A=A 



Define a linear function F by 



n d—1 



df 



a=l i=l dX 



(a) 



A 



(a) 



A=A 



Then Eq. fl7H|) implies that F{Li^i^...i^) must be the same for all [iii2 . . . in] G A^, where 
is the set of n-tuples such that Ci^i^...i^ 7^ for the given state. Let this value be v. Then 
the equation -F(A) = v determines a hyperplane containing all lattice E-points Li^i^...i^ such 
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that [iii2 ■ ■ ■ in] ^ N. This imphes that the convex set is contained in a hyperplane. 
Thus, C Vat with N C N, where Vn is a hyperplane. Hence, by Lemma EHl we find 
A E Vff C. Vn- Thus, all E-points on /(A) =0 are on a hyperplane Vn for some N. Since 
there are only a finite number of such hyperplanes, and since the connected hypersurface 
/(A) = has a tangent hyperplane everywhere, this hypersurface must coincide with a single 
hyperplane Vn- □ 
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